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Abstract 



We consider the problem of 2D supersonic flow onto a solid wedge, or 
equivalently in a concave corner formed by two solid walls. For mild 
corners, there are two possible steady state solutions, one with a strong 
and one with a weak shock emanating from the corner. The weak shock is 
observed in supersonic flights. A long-standing natural conjecture is that 
the strong shock is unstable in some sense. 

We resolve this issue by showing that a sharp wedge will eventually pro- 
duce weak shocks at the tip when accelerated to a supersonic speed. More 
precisely we prove that for upstream state as initial data in the entire do- 
main, the time-dependent solution is self-similar, with a weak shock at the 
tip of the wedge. We construct analytic solutions for self-similar potential 
flow, both isothermal and isentropic with arbitrary 7 > 1. 

In the process of constructing the self-similar solution, we develop a large 
number of theoretical tools for these elliptic regions. These tools allow us 
to establish large-data results rather than a small perturbation. We show 
that the wave pattern persists as long as the weak shock is supersonic- 
supersonic; when this is no longer true, numerics show a physical change 
of behaviour. In addition we obtain rather detailed information about 
the elliptic region, including analyticity as well as bounds for velocity 
components and shock tangents. 

Contents 

1 Introduction 

1.1 Background 

1.2 Numerical results 

1.3 Main result 

1.4 Related work 

1.5 Overview 

1.6 Notation 

2 Potential flow 

2.1 Unsteady potential flow 

2.2 Self-similar potential flow 

2.3 Shock conditions 



2.4 Pseudo-normal shocks 

2.5 Shock polar 

2.6 Shock-parabolic corners with fixed vertical downstream velocity 

Maximum principles 

3.1 Common techniques for extremum principles 

3.2 Density in the interior 

3.3 Velocity components in the interior 

3.4 Velocity components on the wall 

3.5 Velocity components at shocks 

3.6 Pseudo-Mach number at shocks 

3.7 Density at shocks 



Construction of the flow 

Problems 

Approach 

Parameter set 

Function set and iteration 

Regularity and compactness .... 
Pseudo-Mach number control . . . 

Second derivatives on arcs 

Arc control 

Corners moving along arcs 

Corner bounds 

Density bounds and shock strength 
Velocity and shock normal control 

Fixed points 

Leray-Schauder degree 

Existence of fixed points 

Construction of the entire flow . . 



4 


1 


4 


2 


4 


3 


4 


4 


4 


5 


4 


6 


4 


7 


4 


8 


4 


9 


4 


10 


4 


11 


4 


12 


4 


13 


4 


14 


4 


15 


4 


16 



Appendix |95 

5.1 Regularity in 2D corners |95 



5.2 Free boundary transformation llOl 



1 Introduction 

1.1 Background 

Gas flow onto a solid wedge, like forward edges of airplane wings or engine 
inlets, is a fundamental problem for aerodynamics (see Figure[2]). An equivalent 
problem is flow in a convex corner of an otherwise straight wall (see Figure [T]). 
For supersonic flow and sufficiently small t, it is well-known that this problem 
has steady solutions with a straight shock emanating from the corner, separating 
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Figure 1: Two steady solutions of supersonic flow along an infinite wall with a 
corner. 



Figure 2: Reflect either part of Figure[T] across the upstream wall. There are four 
combinations; each is a steady solutions to supersonic flow onto a solid wedge. 
The one with two weak shocks is the one observed in nature and numerics. 

two constant-state regions ( "upstream" and " downstream" ) . The shock is the 
consequence of compression of the gas by the downstream wall. 

A longstanding open and puzzling problem is that, for r close to (correspond- 
ing to a sharp wedge resp. a mild corner), there are two possible steady solutions 
of the corner flow, one with a strong and one with a (comparatively) weak shock 
(see Figure [T]). Both shocks satisfy the entropy conditior0. However, only the 
weak shocks are observed in nature. To quote |CF48j : "The question arises 
which of the two actually occurs. It has frequently been stated that the strong 
one is unstable and that, therefore, only the weak one could occur. A convincing 
proof of this instability has apparently never been given." 

The goal of the present paper is to understand this. 

For many purposes, in particular for many questions concerning flow around 
airplane wings, viscosity, heat conduction and kinetic effects can be neglected. 
It is natural to consider inviscid models, such as the full or isentropic compress- 
ible Euler equations or compressible potential flow. The appropriate boundary 
condition at solid surfaces is the slip condition: the gas velocity is tangential. 

In each model the shock and its upstream and downstream states satisfy the 
Rankine-Hugoniot relations, a system of nonlinear algebraic equations. These 

^ There is a third shock that violates the entropy condition. 
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Figure 3: Curve: downstream velocities for all possible shock normals. The 
left branch corresponds to compression shocks, the right branch to unphysical 
expansion shocks. The black + is exactly sonic, points left of it are subsonic, 
right of it supersonic. Solid line: there are three solutions; strong shock K, 
weak shock L, unphysical. Dotted line: still three solutions, but the weak shock 
is supersonic-subsonic now. Dashed line: critical angle; weak and strong shock 
coincide. Dashed-dotted line: t above critical, no physical solution. 

relations determine the shock polar: the curve of downstream velocities that 
result when varying the shock normal while holding the shock steady and the 
upstream velocity and density fixed (see Figure [3]). 

In Figure [3] the upstream velocity is (1,0), labelled /. Possible downstream 
velocities are intersections of the shock polar with the ray at a counterclockwise 
angle r from the positive horizontal axis. Obviously for small t there are three 
intersection points. The leftmost intersection, called K, corresponds to the 
strong shock. The rightmost is an unphysical expansion shock which need not be 
considered. The middle point, called L, is the weak shock. (The shock normals 
are parallel to the difference between downstream velocity (intersection point) 
and upstream velocity (1,0).) 

For r 1 0, the strong shock approaches a normal shock, whereas the weak shock 
vanishes {L approaches /). 

There is a critical angle t = where L and K coincide; for larger r no steady 
entropy-satisfying shock can be attached to the wedge tip resp. corner. 



4 



The black + on the shock polar indicates a downstream state that is exactly 
sonic (Macli number AI = 1). Polar points left of it are subsonic {M < 1), 
polar points right of it arc supersonic. In numerical experiments, the weak 
shock detaches from the corner/wedge tip when its downstream changes from 
supersonic to subsonic (e.g. by increasing r). 

In |EL05bl |EL06| we have reported on numerical experiments: to our surprise 
and somewhat contrary to the aforementioned conjecture, both the strong and 
weak shocks arc time-asymptotically stable under large, compactly supported 
perturbations. Instead, the strong shock is unstable under (generic) perturba- 
tions of the downstream state at infinity; depending on the perturbation cither 
the weak shock appears or the shock detaches from the wedge tip/corner entirely. 
It may be possible to obtain a strong shock in very special cases, for example 
by placing a perfectly feedback-controlled nozzle somewhere downstream. 

The weak shock is stable under both kinds of perturbation. 

While various conjectures and empirical observations have been made regarding 
weak vs. strong shock, previously no mathematical arguments for either were 
known. To obtain one, we devise an "unbiased" test: at time t = 0, fill the 
entire domain with upstream data; check which shock appears for t > 0. 

In numerics, the weak shock appears spontaneously (see Figure 2]). Motivated 
by this, we construct an analytical solution. 

An equivalent experiment is to accelerate a solid wedge in motionless air instan- 
taneously to supersonic speed. More generally, if a finite wedge is accelerated 
from rest at time to a fixed supersonic speed at time e <^ L/c {L wedge 
length), we may expect the solution to be a good approximation for times t in 
the scale e -^t <^ L/c. 

1.2 Numerical results 

Figure S] shows the flow pattern (density) solving our test problem, for some 
positive time. Here 7 = 7/5, M/ = 2.94 and t = 10°. vj is horizontal from left 
to right. Blue is pi] green, yellow and red are successively larger densities. 

A straight shock (blue to red) emanates from the wedge tip. Calculation shows 
that it is the weak shock. There is another straight shock on the right (blue 
to orange), parallel to the downstream wall. Below each straight shock lies a 
constant region. Both shocks are connected by a curved shock, with a nontrivial 
(elliptic) region below. 

The flow pattern is self-similar: density and velocity are constant along rays 
X = £,t,y = i]t for fixed ^, rj. This can be visualized as t being a "zoom" 
parameter, with t = corresponding to "infinitely far away" and i t 00 to 
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Figure 4: Left: numerical solution of the wedge flow problem. Right: Structure. 
Three trivial hyperbolic regions I, L, R are separated by straight shocks. They 
enclose a nontrivial elliptic region bounded by a curved shock and two parabolic 
circle arcs centered in resp. vn with radius Cl resp. cr. Density and velocity 
are functions of ^ = x/t. 77 = y/t only. 

"infinitely close to the origin" (wedge tip resp. wall corner). In particular the 
flow structure is the same for all times. 

Figure [S] shows the elliptic region in more detail. In the top picture density is 
shown; its minimum in the elliptic region is attained at the shock (Proposition 
13.7.11 will show that the minimum is a "pseudo- normal" point). In the middle 
picture the velocity tangential to the downstream wall is shown; the bottom 
picture shows velocity normal to the downstream wall. 

The diamond in the center of the bottom domain boundaries indicates the origin 
in self-similar coordinates, where we use standard coordinates (Figure [12]) ■ 

It should be emphasized that numerical computations only suggest the structure 
of the solution. For instance, it is not clear that the constant states L and R 
extend to the pseudo-sonic circles Pl and Pr. Although in one dimension with 
viscosity, some techniques can convert a numerical solution with sufficiently 
small residual into an existence proof for an exact solution (see e.g. }JY98| 1. 
only partial results are available in multiple dimensions without viscosity (see 
|Kuz75[ ILW60[ IE1107| ) ; it is not clear whether a general result is even true (see 
pl06] l. 

1.3 Main result 

To obtain a mathematical argument, we construct the self-similar solution ex- 
actly rather than numerically. We use compressible potential flow as model: 

Theorem 1. Let r e (0, f), Mi,pi,ci G (0, 00), 7 e [l,oo); set vj := 
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Figure 5: Elliptic region. From top to bottom: p, v tangential, and v normal 
to the wedge surface. Corner/wedge tip is intersection of left shock and bottom 
domain boundary (left outside diagrams). Diamond indicates origin in Figure 
[T2] coordinates; velocities are relative to it. 
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(M/C7,0). Define the wedge 

W := {(x^y) e IR^ : y < xtanr} 

and 

n (0,oo) X CW. 
Assume the following conditions are satisfied: 

1. Unsteady potential flow with ^-law pressure admits a steady straight shock 
with upstream data pj and vj and downstream velocity vl and sound speed 
Cl so tha^ 

1{vi,vl) = T. 

2. The shock is supersonic- supersonic: 

Ml ^ > 1. 

CL 

3. Of the two intersection points of the shock with the circle dBc^{vL), let 
^£ be the one closer to the corner (origin). Let the R shock be the unique 
shock parallel to vl, with upstream data pi and vj, downstream sound 
speed cr and downstream velocity vr parallel to vl as well. Of the two in- 
tersection points with dBcp^{vR), let ^'^ be the one farther from the corner 
(see Figure\^. We require that 

{Line segment from to H Bci{vj) = (b. (1.3.1) 

Then there is a weak solution (see Remark ] L 3.1]) (f> G C^'^{Q?) of 

unsteady potential flow in VL, (1.3.2) 
V0 • n = on dW, (1.3.3) 

p = V</) = w/ fort = (1.3.4) 

In addition to existence, detailed results about the structure of the weak solution 
can be obtained (see Remark 14.16. 2p . At this point we emphasize only that 
each solution consists, in some neighbourhood of the origin, of the weak shock 
separating two constant-state regions. 

Remark 1.3.1. See Scction r2.1l for introduction and precise definition of potential 
flow. By weak solution wc mean that 

V0(O, x) = VI for a.c. xeCW (1.3.5) 

is the counterclockwise angle from first to second vector, ranging from to 2n. 
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Figure 6: Solutions are constructed for all cases that satisfy the condition p.3.ip : 
the dashed circle and line must be separated. A shock has positive strength if 
and only if it does not touch the circle. 

and 

/ p-dt + pV(j) -V-d dx dt+ / z?(0, x, y)pi dx ^ 
Jn Jew 

for all test functions i9 £ C^m- 

(For e C°'^(ri), the velocity V0 is a.e. well-defined on {0} x ZW, but (pt and 
hence p may not be well-defined.) 

Remark 1.3.2. As Remark 14.16.11 shows, there is a large set of tip shocks and 
parameters that satisfy the conditions of Theorem [1] 

The first and second condition are physically necessary, not technical limitations. 
If the first is violated (for Mj < 1 or large r), there is no straight steady physical 
shock attached to the corner at all. If either of them is violated, numerical 
experiments show a flow pattern with a shock detached from the corner, moving 
upstream (left). 

The third condition is technical. It is needed in some cases to prove the shock 
docs not vanish (which is never observed in numerics); none of the other esti- 
mates requires this condition. We expect that the condition will be removed 
with some additional analysis. 

It should be emphasized that the theorem and its proof are global in nature: 
large parameter changes are possible. 

Remark 1.3.3. Incidentally we also solve the problem for asymmetric wedges, 
as long as both sides allow a supersonic-supersonic weak shock and as long as 
(|1.3.1[) is satisfied on both sides. 

1.4 Related work 

|CF48[ Section 117, 122 and 123] explain in detail shock polars and the corner 
flow problem. Despite its age, its discussion of weak vs. strong shock is still a 
good reflection of the state of prior research. |FT68j is another useful reference. 
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CanicKK02| consider the classical problem of regular reflection of a shock by 
a symmetric wedge; this problem, like ours, has a self-similar solution. They 
consider the unsteady transonic small disturbance equation as model. [Zhe06j 
studies the same problem for the pressure-gradient system. The monographs 
|Zhe01[ ILZY98j compute various self-similar flows numerically and present some 
analysis and simplifled models. 

|CF03| prove existence of small perturbations of a plane shock in steady potential 
flow. 

|Che03j constructs steady solutions for 3D cones rather than 2D wedges. jLL99) 
discuss stability of 3D flow past a perturbed cone; |CL05| show existence and 
linear stability in the case of the isentropic Euler equations. |CZZ06| study 
existence and stability of supersonic flows onto perturbed wedges with attached 
shocks; the introduction gives a detailed discussion of previous work. 

So far the only other paper that proves global existence of some nontrivial time- 
dependent solution of potential flow is |CFj : they construct exact solutions for 
regular reflection, assuming sufhciently sharp wedges. 

1.5 Overview 

In Section [2] we give an introduction to unsteady potential flow. We derive self- 
similar potential flow, discuss its shock conditions and analyze the properties of 
shocks in detail. 

In Section [3] wc discuss a collection of maximum principles for elliptic regions of 
self-similar potential flow. Some of these identify circumstances in which certain 
quantities (density, ...) can or cannot have maxima or minima in the interior. 
Other results discuss local extrema at solid (slip condition) walls and flnally 
shocks with a constant-state hyperbolic region on the other side. 

Since the hyperbolic regions are trivial (see Figured]), the heart of the problem 
is the construction of the elliptic region. This is accomplished in Section 21 
Readers interested in more overview should go to Section [121 where all proof 
steps are surveyed. 

A crucial ingredient are the maximum principles from Section [3l combined with 
ODE-type arguments at the parabolic arcs in Sections 14.61 to 14.101 a-nd tech- 
niques to control shock location and normals (Sections 14. Ill and I4.12p . Section 
14.161 combines the elliptic region with its hyperbolic counterparts to construct 
the full flow pattern. The remaining sections are standard but delicate applica- 
tions of nonlinear elliptic theory. Some literature results, such as regularity in 
corners and at free boundaries, need extension which is done in the Appendix. 
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1.6 Notation 



For the most part we use standard notation. Subseripts and superscripts may 
denote tensor indices, partial derivatives or powers, depending on the context. 

y) is the counterclockwise angle from x to y. For x = (a;i,X2), y = (2/1,2/2), 

x-^ := {-X2,xi) 

{counterclockwise rotation by 90°), 

xxy:^ xiy2 - X2yi- 

x^ is the rank-one matrix xaf^ whereas |xp is the norm. Correspondingly, 
= VV"'" is the Hessian {not the Laplacian). 

Normals n are outer normals to a domain, except on the shock S (defined later) 
where they are downstream, so usually inner. Tangents t are always defined as 
t := n-L. 



2 Potential flow 

2.1 Unsteady potential flow 

We consider the isentropic Euler equations of compressible gas dynamics in d 
space dimensions: 

pt + \/-{pv)^0 (2.1.1) 

d 

{pv)t + Y.^pv'v),. + V{p{p)) = 0, (2.1.2) 

Hereafter, V denotes the gradient with respect either to the space coordinates 
X = (x^,a;^, • • • ,x'^) or the similarity coordinates t~^x. v ~ (w^,u^, • • • ,1)'^) is 
the velocity of the gas, p the density. p{p) pressure. In this article we consider 
only polytropic pressure laws (7-laws) with 7 > 1: 

Pip)-'^(^Y (2.1.3) 
7 \PoJ 

(here cq is the sound speed at density po). Many subsequent results extend 
with little or no change to 7 < 1 or to general pressure laws, but in special cases 
some steps require more work or break down entirely. To keep the presentation 
simple we don't strive for generality with respect to pressure laws. 

For smooth solutions, substituting (|2.1.ip into (|2.1.2p yields the simpler form 

Vt+v-\/'^v + \/{Tr{p)) = 0. (2.1.4) 
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,-1 7>1 



Here tt is defined as 

ilog(p/po), 7=1- 
This TT is C°° in p G (0, oo) and 7 G [1, cxj) and has the property 

Pp 



If wc assume irrotationality 



v] = vl 



(where z, = 1, . . . , d), then the Euler equations are reduced to potential flow: 

V = V3(j) 

for some scalar potential function (j). For smooth flows, substituting this into 
(|2.1.4p yields, for z 1, . . . , d, 

= cj),t + • Vc/. + ^(p), = (0t + + ^{p))^. 

Thus, for some constant A, 

p = .-i(A-0,-^). (2.1.5) 

Substituting this into (|2.1.ip yields a single second-order quasilinear hyperbolic 
equation, the potential flow equation, for a scalar field 0: 

{p[c^u |V(/.|))^ + V • {p[c^u |V0|)V0) = 0. (2.1.6) 

Henceforth we omit the arguments of p. Moreover we eliminate A with the 
substitution 

^ ^ 0, (j){t, x) ^ (j){t, x) - tA 
(so that (pt (pt — A). Hence we use 

p = ^-i(-0t-i|V0|2) (2.1.7) 

from now on. 

Pp 

(--^)' = (-pr-(^)-^ = ^ (2.1.8) 



^We consider simply connected domains; otherwise (j) might be multivalued. 
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the equation can also be written in nondivergence form: 

d 

(jjtt + 2\7(j)t • V(/) + ^ - c^Acj) = (2.1.9) 

(|2.1.9p is hyperbolic (as long as c > 0). For polytropic pressure law the local 
sound speed c is given by 

c' = c^ + (7-l)Mt-^|V</>|^). (2.1.10) 
2.2 Self-similar potential flow 

Our initial data is self-similar: it is constant along rays emanating from x = 
(0,0). Our domain CVF is self-similar too: it is a union of rays emanating from 
{t,x,y) = (0,0,0). In any such situation it is expected — and confirmed by 
numerical results — that the solution is self-similar as well, i.e. that p, v are 
constant along rays x = emanating from the origin. Self-similarity corre- 
sponds to the ansatz 

0(t,f) := tV(f), ^■.= t-^x. (2.2.1) 
Clearly, € C°'^{n) if and only if ?A G C°'^{ZW). This choice yields 
v(t,x) = V(/)(t,f) = Vipit-'^x), 

p{t,x) = 7T-\-cf,t - i|v0n = n-\-iJ+f w - ^ivv-n. 

The expression for p can be made more pleasant (and independent of .f) by using 

x(e):=VXO-^l?T; 

this yields 

P = ^-\-X-INx\')- (2.2.2) 

Vx = Vt/i — ^ is called pseudo-velocity. 
(|2.1.6p then reduces to 

V-(pVx) + 2p = (2.2.3) 

(or +dp, in d dimensions) which holds in a distributional sense. For smooth 
solutions wc obtain the non-divergence form 

{c'l - VxVx^) : V^x = (c2 - xDxii - ^XiXvXiv + (c' " x')x,, = |VxP - 2c' 

(2.2.4) 
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Another convenient form is 



(2.2.5) 

Here, (|2.1.10p for polytropic pressure law yields 

^=cl + {^-l){-X-\\Vx?) (2.2.6) 
Remark 2.2.1. (|2.2.3p inherits a number of symmetries from (|2.1.ip . (|2.1.2I) : 

1. It is invariant under rotation. 

2. It is invariant under reflection. 

3. It is invariant under translation in ^, which is not as trivial as translation in 
x: it corresponds to the Galilean transformation v ^ v -\- vq, x <— x — vot 
(with constant vq G IR'') in {t,x) coordinates. This is sometimes called 
change of inertial frame. 



(|2.2.4[) is a PDE of mixed type. The type is determined by the (local) pseudo- 
Mach number 

L:=l^, (2.2.7) 
c 

with < L < 1 for elliptic (pseudo-subsonic), L = \ for parabolic (pseudo- 
sonic), L > 1 for hyperbolic (pseudo-supersonic) regions. 

The pseudo-Mach number L can be interpreted in a way analogous to the Mach 
number M: consider a steady solution of the unsteady potential flow equation. 
Loosely speaking, in an Af < 1 (subsonic) region a small localized disturbance 
will be propagated in all directions, whereas in an M > 1 region it is propagated 
only in the Mach cone. L < I and L > 1 are analogous, except that we study 
the propagation of disturbances in the unsteady potential flow equation written 
in (t^x/t) coordinates, rather than {t,x). 

There is no strong relation between M < 1 and L < 1: consider two constant- 
state (hence steady and selfsimilar) flows with zero velocity {M = constant) 
resp. supersonic velocity v {AI > 1 constant). Each flow has L = in the point 
= resp. ^ = V and L t c« as ^ t c», so there are examples for each of the 
four cases M, L < 1, M < 1 < L, L < 1< AI and 1 < L, M. 

While velocity v is motion relative to space coordinates x, pseudo- velocity 

f := Vx 

is motion relative to similarity coordinates ^ at time t = 1. 
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The simplest class of solutions of (|2.2.4p are the constant- state solutions: ip 
afhne in ^, hence v, p and c constant. They are elliptic in a circle centered in 
^ = V with radius c, parabolic on the boundary of that circle and hyperbolic 
outside. 

Convention 2.2.2. If we study a function called (e.g.) x, then 'ip, p, L etc. 
will refer to the quantities computed from it as -ip, p, L are computed from x 
(e.g. V' = X + ilCP)- We will tacitly use this notation from now on. 



2.3 Shock conditions 

Consider a ball U and a simple smooth curve S so that U = ?7" U SUlf^ where 
U'^jU'^ are open, connected, and S,U'^,U'^ disjoint. Consider % : f/ — > IR so 
that X = X"''' in U"^'^ where x"'"* e C^{U^). 

X is a weak solution of (j2.2.3p if and only if it is a strong solution in each point 
of U- and U+ and if it satisfies the following conditions in each point of S: 

X"=X', (2.3.1) 
n • (p"Vx" - p'^Vx'^) = (2.3.2) 

Here n is a normal to S. 

(|2.3.1[) and (|2.3.2p are the Rankine-Hugoniot conditions for self-similar potential 
flow shocks. They do not depend on ^ or on the shock speed explicitly; these 
quantities are hidden by the use of x rather than ip. The Rankine-Hugoniot 
conditions are derived in the same way as those for the full Euler equations (see 
|Eva98( Section 3.4.1]). 

Note that (|2.3.ip is equivalent to 

V'" = ip'^. (2.3.3) 
Taking the tangential derivative of (|2.3.1[) resp. (|2.3.3[) yields 

9x" dx'^ 

or dyj" 

The shock relations imply that the tangential velocity is continuous across 
shocks. 

Define (zj, z^) :— z„ :— Vx" and (uj, v!^) :— Vu VV'". Abbreviate := Zu -t, 
zl^ := Zu ■ fi, and same for v instead of z. Same definitions for d instead of u. 
We can restate the shock relations as 

Puz: = pdz2. (2.3.6) 

4 = 4- (2.3.7) 



(2.3.4) 
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Using the last relation, we often write z* without distinction. 

The shock speed is a = ^ ■ n, where ^ is any point on the shock. A shock is 
steady in a point if its tangent passes through the origin. We can restate (|2.3.6p 
as 

P«< - PdVa = ct(Pu - Pd) 
which is a more familiar form. 

We focus on p^, Pd > from now on, which will be the case in all circumstances. 
If p„ = in a point, we say the shock vanishes; in this case = z" in that 
point, by (|2.3.7|) . In all other cases z^,z^ must have equal sign by (|2.3.7p : we 
fix n so that z2,z^ > 0. This means the normal points downstream. The shock 
is admissible if and only if p„ < pd which is equivalent to z^ > z^. 

A shock is called pseudo-normal in a point ^ if z* = there. For ^ = 0, this 
means that the shock is normal (w' = 0), but for 7^ normal and pseudo- 
normal are not always equivalent. 

It is good to keep in mind that for a straight shock, pd and Vd are constant if 
Pu and Hu are. Obviously Zd may vary in this case. 

Remark 2.3.1. The Rankine-Hugoniot conditions for the original isentropic Eu- 
ler equations cannot be used for potential flow: even if the flow on one side of a 
shock is irrotational, the flow on the other side has nonzero vorticity for curved 
shocks. 

2.4 Pseudo-normal shocks 

Here we study the consequences and solutions of the shock relations in a pseudo- 
normal point. We state all results for pseudo-velocities z and for pseudo-Mach 
numbers L because moving shocks are ubiquitous in this article. For better 
intuition the reader may bear in mind that z = v and L = M if the shock 
is steady, i.e. passes through the origin. In fact by Remark 12.2.11 weak and 
entropy solutions of self-similar potential flow are invariant under translation 
and rotation, so we may always consider translating the shock so that it becomes 
steady, which does not change z, p, L whereas v is changed only by a constant 
vector. Hence the behaviour of arbitrary shocks is entirely determined by those 
of steady shocks. 

L" := z"/c and := z*/c will be referred to as normal resp. tangential pseudo- 
Mach number. (|2.3.ip and (|2.2.2p imply 

P. = --^Wp„) + ^~^). (2.4.1) 
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It is apparent that (|2.4.ip reduces to 

P, = n-\n{p^) + 4^ - 4^). (2.4.2) 

Combined with (|2.3.7p there is a direct relation between normal velocities, in- 
dependent of the tangential velocities. 

For polytropic pressure laws we may use a rather convenient simplification: 
there is an explicit relation connecting L", LJJ, independent of p^, z^. 

Lemma 2.4.1. For > 0, (§4^ and 12.3.6]} are equivalent to 

g{Ll)^g{L-,), (2.4.3) 

„2 , 2 ^ 



Note that 



X — 21oga;, 7=1. 



|^ = ^(x-x-i)x-2^, (2.4.5) 
ax 7+1 

d^g 4 



^-(:^((3-7) + (37-l).->--^. 



Moreover 



as well as 



r f TJ] \ ~ 



(2.4.6) 



gix) /^"^ ^ > I asxlO, (2.4.7) 

[-logx, 7=1 

4 

17(0;) ^ xT+i as x 1 00, (2.4.8) 



(2.4.9) 



(2.4.10) 



Proof. (|2.3.6p can be written 

PuL'^Cu = PdLaCd 
=> {p^)('>+'y^L: = p[^+'^/'L-,; (2.4.11) 

this yields (|2.4.10p which yields (PX^)) . 



17 



Consider 7 > 1, so that = Cg + (7 — l)7r(p). (|2.4.2p can be transformed to 

2 7-1 2 7-1 

Substitute ()2.4.9p to obtain (|2.4.3p . 
For 7=1, (PX^ is 



Po 2 " ' Po 

P« 2 

using (|2.4.10p we obtain (|2.4.3p . □ 

Proposition 2.4.2. There is an analytic strictly decreasing function = 
^di^u)' which is its own inverse, so that the shock relation {2.4-^ is solved for 
all L" S (0,cx3). For L" 7^ 1, t/ie onZj/ other solution of \2.4-.S^ is the trivial 
one: Ui = L" . For LI] — 1, both coincide. 



r 00, 0, 

= 1, K = 1, 

-{K)^\ ^ K^^. 7>i, 

~exp(-(^), Too, 7=1. 



The resulting shock is admissible if and only if > 1 . 



(2.4.12) 



dLV, 



93 

dLl] 



|iS=iS(=i) 



Ll -l/Ll (L 



rn 



= -1, 



4' 



(2.4.13) 
(2.4.14) 
(2.4.15) 



For L" > 1 > L^, we have pd > Pu o,nd < 2". Pd,Cd are strictly increasing 
in L" for p„ fixed. For 7 > 1, > c„ as well, and Cd is strictly increasing in 
Lu for Pu fixed. 



Proof. Assume > 1. From (|2.4.5p it is obvious that |f (a;) > for x > 1. 
Hence g(L^) > 5(1); moreover (|2.4.3p cannot have more than one solution LJJ 
in [l,c»] for fixed > 1; in fact L'2 = L" is the unique solution. For a; < 1, 
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< 0, so (|2.4.3p cannot have more than one solution L'2 in (0, 1). It must 
have one, though, because .9(0+) = +oo > g{L") > g{l). 

For iJJ < 1, the existence of a nontrivial solution S (I7O0) is obtained from 
the previous case by analogous arguments. For L" = 1, the sign of dg/ dx rules 
out any other solutions. Since (|2.4.3p is symmetric in L",L'2, it is obvious that 



is its own inverse. 



The trivial solution branch = i" is obviously smooth; by the implicit func- 
tion theorem the the nontrivial branch is analytic away from = 1. In L" = 1 
there is a degeneracy which has to be analyzed by inspecting the Hessian of 

hiL^,L",) := g{Ll) - g{L-^) = 0: 



A 



WW 



WW 



d: 



9 (Tn\ 









1 12.4.61 



8 
7 + 1 







-8 
7+1 



A is an invertible indefinite matrix; the solutions of nF Aw = are w ~ (1, 1) 
and w = (1,-1). The classical Morse lemma (see e.g. }Smo94| Lemma 12.19]) 
shows that in a small neighbourhood of (LJ^, LJJ) = (Ij l)i the solution of (|2.4.3p 
form two analytic curves that intersect in (1,1) with tangents (1,1) (trivial 
branch) and (1,-1) (nontrivial branch). The latter yields (|2.4.14p . 

The remainder of (j2.4.12p follows from the asymptotics of .g (see p.4.7|) . (|2.4.8|) ). 



From (|2.4.10p and (|2.4.13|) we see that pd > Pu for < 1 < and that pa is 
an increasing function of L" for pu held fixed. Clearly the same applies to if 
7 > 1, and (because of (|2.3.6p ) < z". This means the shock is admissible for 
> 1. □ 

Remark 2.4.3. In the remaining arguments we always assume that the shock is 
admissible (or vanishing) and ignore the branch = U^- 



Proposition 2.4.4. If we hold c^, Pu fixed: 



In particular 



Therefore 



dz^ 

dzl: 



dzl 



< 



7-1 
7+1 



< 1. 



< 



7-1 
7+1' 



(2.4.16) 



(2.4.17) 



I£L_1 < J 1 

Cu 7+1 



(2.4.18) 



Remark 2.4.5. (|2.4.17p is not very tight (we can show < for 7 < 3), but 
sufficient for our purposes. 
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dcd _ d fUl\^^ cu fL^\^' i~l f^ LldU^ 



Proof. We use ([^X^ : 

dcd d 



so 



dLl 7+1 

zJJSzJ pafL^\^f 2 , 7-1 



Then 



z^^9z» p„ V^S/ V7 + 1 L2 dLl 7 + 1 
12X101 2 dL^ 7-1 7-1 

7 + 1 L;J 7+1 7+1' 

>0 >0 <0 

Integrating (|2.4.17p from — z"!^ — Cu for a vanishing shock, we obtain (|2.4.18p . 

□ 

Proposition 2.4.6. Consider a shock with velocity a := ^ ■ ft. Our convention 
z" > requires a < u". Vary a while holding n and Vu fixed. Then: 

?^>l-^>0; (2.4.19) 

^ < 0. (2.4.20) 
oa 

Proof. For moving shocks, v^J: = z^ + (J and z^ = — cr, so 



> 1 - = > 



da da dzj^ 7 + 1 7 + 1 

and 

dpd ^ _dpd ^ Q 
da dz^ dL^^" fTim 

2.5 Shock polar 

Here we prove only the resuhs needed for our purposes. 



□ 
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Proposition 2.5.1. Consider a fixed point on a shock with upstream density 
Pu and pseudo-velocity Zu held fixed while we vary the normal. Define (3 := 
Z(z„,n). pd is strictly decreasing in |/3|, whereas Ld, \zd\ are strictly increasing. 
Cd is strictly decreasing for ^ > 1, constant otherwise. Moreover 

{dpvd) ■ n = (dpzd) . r? = z*(^ - l), (2.5.1) 
(dpvd) ■ t = (dpzd) ■ t = 4 - z^. (2.5.2) 

If = (-^uiO) 'with z'^ > 0, then z'^ is increasing in |/3|. 
Proof. 

dfiZd = dpiz^ri + /t) = ^dpz^n + z2t- z'n + dpzH 
oz^ 

\ OZ ^^ / 

This is (PXT|) . (PX^ . using that ^ is fixed. For /3 > 0, z* = -|Zi,|siii/3 is 
strictly decreasing (and negative). 

\dp{\zd?) = \dp{{z-,f + izY) = ^2^dpz- + z'dpz' 



d 



dz^ 

<(7-l)/(7+l) 

by (|2.4.17p . so \zd\ is strictly increasing. Then by (|2.4.ip 
Pd = ^"'(7r(p„) + i(|i;p_|^^|2)-)^ 

is strictly decreasing, as is Cd (except constant for 7 = 1), so Ld = \zd\/cd is 
increasing. 

For /3 > 0, ([2X3)) yields 

dp{zi) - - 1)^ - (^5_-^£)sin^ > 

<o ^° 

by (12.4.171) . 

/3 < is analogous by symmetry. □ 
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Figure 7: L values downstream of a shock. 



2.6 Shock-parabolic corners with fixed vertical downstream 
velocity 

The corners of our eUiptic region (see Figure [T^ are points on shocks where 
Ld = 1 (or Ld = Vl ~ £ if regularized); on the other hand = (if the x 
direction is tangential to the wall). We study such shocks in detail. 

Proposition 2.6.1. Consider a straight shock (see Figure^ passing through a 
point ^. The shock is pseudo-normal in the point S,m = ^+t-(vd—£,)'t cind pseudo- 
oblique in every other point. S,m is the closest point on the shock both to Vu and 
to Vd- The circle with center Vu and radius Cu does not intersect the shock. 
The downstream flow has L < VI — e inside the circle with radius Cdy/i — e and 
center Vd, L = \J\ — e on the circle and L > y/\ — e outside. If < ^\ — e, 
then the circle intersects the shock in the two points £,m ± Cd-\/l — e — [L'^)'^ ■ t. 

Proof. Straightforward to check. ■ t — Vd ■ t ~ Vu ■ t, so \vu ~ £,m\ = \{vu — 
I'm) • n\ = \z:\ > Cu. □ 

Proposition 2.6.2. Consider a straight shock with = 0, v]^ < and down- 
stream normal ft = (sin /3, — cos/3) through ^ ~ (0, 77) (see Figure\^. For every 
(3 G §) there is a unique r] = tjq € \R so that = 0. rjg and the corre- 

sponding downstream data are analytic functions of (3. tjq is strictly increasing 
in \(3\. 

For the shock passing through (0,77g), let ^£ and be the two points with 
Ld = \/l — e, as given by Provosition \2.6.1\ These points are analytic functions 
of p. L^, pd and z" are increasing function^ of (3; and are decreasing 
functions of j3. For (3 £ [0, ^), r/2 i^ strictly decreasing function of (3 with 
range (j^^i ^^o]; where tJq is the tJq for f3 ~ 0, and rj*^^ is some negative constant. 

''AH of these are independent of the location along the (straight) shock. 
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Proof. First regard everything as a function of /S and 77, with Pu held fixed 
and the shock held as passing through (0, 77). In (0, 77): 

= Zun"" + = {< - ^) sin /3 - {vl - 77) cos /3 ^ {ij - vl) cos /3, 
=0 =0 

(2.6.1) 

d^z^ = cos/3, 

«^ = '^^l + (^'rf' - <) - ^'il + K - vl)ny =vl + {z2 - z:)ny = vl + {z^ - z^^) cos /3, 

(97™ (|2 4.171 O 

dv^l = (1 - cos/3 = (1 - ^) cos2 p cos2 /3 > 

^ ^ dzZ^ dzl' 7+1 



(2.6.2) 

tJ^ is an increasing function of 77, so for fixed (3 there can be at most one rj with 
t;^ = 0. For rj ~ v\ ^Cul cos/3 we have L" = 1 in the point (0, 77), so i;^ = < 
there; on the other hand (j2.6.2|) has uniformly lower-bounded right-hand side 
(for every fixed /3e(— f,-|),S0 7;^t +00 if we take r\ \ -l-oo. Therefore there is 
exactly one solution 7/ = 77Q for each /3. 

Now consider /3 G (0, J) first (so z*' — (w^ — 77) sin /3 < 0); the case /3 < is 
symmetric. 

= (1 - 7^)5/3< cos/3 - (z^ - z,") sin/3 = (1 - -^)z* cos/3 - (z^^ - z^') sin/3 
cZjj cz„ 

1 „n n 

^^-v ' >0 V ' >0 

>o >o 
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For the remainder of the proof, fix ij = tJq and consider everything a function 
of p. Consider f3 £ (0, f ) increasing (the other case is symmetric). Then 

2" = i^u 



=z* =— cos;3 

Therefore LJ^ and are also strictly increasing, whereas is strictly decreas- 



■71^ 11^ - 



{Vd - vl)n^ " ( vl - vDn^lny = tan/3 < 



=0 =0 <o 



V 



cos"^ p 

Obviously is strictly decreasing. 

mi ^vl~ zlny = vl + zl cos /? 
i^^g??*/ = sin/3 + dpzl cos/3 = ( - + ) sin/3 

^fWZ^|_ii,^^^™0 ' " (2.6.3) 

So 77M is strictly decreasing. 



= mi - Cd sin/3y^l -e - (L2)2 
Cd is increasing and strictly decreasing, so 77^ is strictly decreasing. Moreover 

Vl= vl +z'^cosP~CdSinP^l - e - {U^f = (i^ cos/3 - sin/3y^l - e - iUi)^)cd 
=0 

L'2 is < 1 — e for P = Q and decreasing in /3 > 0, so it is uniformly bounded above 
away from 1 — e as /3 t f ■ For /3 large enough the rjl expression is negative. 
Thus rjl covers an interval (77* ,775] for some 77^ < 0. □ 



3 Maximum principles 

In this section we derive many a priori estimates for smooth elliptic regions and 
for smooth shocks separating elliptic and constant-state hyperbolic regions. 
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3.1 Common techniques for extremum principles 



Lemma 3.1.1. Let m > 1. If clq, . . . , a„i define a positive semidefinite tensor, 
i.e. 

j=0 ^ 

and if 

Aak + Bak+i+Cak+2=0 Vfc e {0, . . . , m - 2} (3.1.1) 
with constants A, B,C CI IR so that AAC > , then ao = ■ ■ ■ = a„i = 0. 

Proof. AAC > B^ means that A + Bz + Cz^ = has two roots z,z e C - IR. 
Then 

flfe = ^{az'') (fc = 0, . . . ,to) 
where a G I is some Unear combination of ao, ai, so 



k=0 ^ ^ k=0 ^ ^ 



(aizC + vD 



for all S, € IR^. Since 3z 7^ 0, (C,?7) i-^ {z^ + ?/)™ is onto C. The inequality 
cannot be true unless a = 0, so ao ~ ■ ■ ■ — a,„ = 0. □ 

Lemma 3.1.2. Let m > 3 and b = (61,62) e IR^ - {0}, n e IR^ - {0}. // 
ao, . . . , Om € IR satisfy 

m-l ^ s 

E ( ""i. ) (^i"'^ + &2afc+i)$S""'"'' > Ve e IR^ C • ?i > (3.1.2) 

2 



as well as i3.Ll\) with constants A,B,Ce IR so that AAC > B^ , then 



flfc^O Vfc e {0,...,m}. 

Proof. AAC > B'^ means that A + Bz + Cz^ = has two roots z,z e C - IR. 
The general solution of ()3.1.ip is 

ak='Si{az'^) (fc = 0, ...,m). 

where a € C is some linear combination of ao and ai. Substitute this into 
((3X21) : 

771— 1 ^ 

\\ rn — l — k 



fc=0 

ni—l 



k ^k^m—l — k 



k=0 

5R(a(6i+62z)(ze + ?7)'""') VCelR^, e-^>0 (3.1.3) 
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We may use £, ■ n > by continuity, which defines a closed halfplane of IR^. 
2: G C — IR, so the range of 77) i-^- ^ + Z77 is a closed halfplane of C. The range of 
(^,77) 1-^ {£, + ZT])"^~^ is all of C because m— 1 > 2. Moreover 61+622; 7^ because 
either 62 = 0, then h + 632 = 61 7^ 0, or 62 ^ 0, then 9(6i + 632) = b2Q{z) ^ 0. 
Thus the range of (^,77) i-^- a(6i + 62z)(2;^ + 77)™"^ is all of C, contradicting 
p.l.3p . unless a = which means = for all fc = 0, . . . , m. □ 

Lemma 3.1.3. Consider an open set U and a point ^0 G U ■ Assume that there 
is a n ^ (quasi an inner normal) so that for every ^ with ft ■ ^ > there is a 
6 > with 

Uo + te:t e (0,<5)} C U. 

Let -0 be an analytic solution of i2. 2. 5\) in U , so that L < 1, p > and D'^ip = 
in ^0- Then 

a 

(with f a C°° function of its arguments) cannot have an extremum in ^q, unless 
ip is linear or 

'|^Vx = mio- (3.1.4) 



d(yijj) da 

Proof. We use dot notation (= etc.) to indicate relations holding only in ^q- 
We show by complete induction over /c = 3, 4, . . . that D'^i/j = 0. Induction step 
(3, . . . , fc — 1 ^ fc > 3): for 7 = 0, . . . , fc — 2 take 9f Sj "^""* of the equation. This 
yields 

(c'/ - Vx') : y^did^^-'^ = (3.1.5) 

because all other terms contain at least one component of D'^i/;, . . . , D''~^ip as 
factor, hence vanish. 



We may exploit that D^x = D^i^ for j > 3. 

2 



2 2 



i=l 2 

and for multiindiccs a with 3 < |a| < fc, 

5"(x+2lvxn = E^^xa,9"v, 

so for all 2 < |a| < fc 

9"(x + ^|Vxn • Vx-V9"7/-, 
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Thus for j = 0, . . . 1: 

+ terms with D'^ip, . . . , D'^^^ijj components as factor 

because D'^tp, .. ., D'^~^ip = 0. 

In a similar way we obtain D^{f ) ^ for ^ ~ 2, . . . , /c — 2, and D{f) = already 
by assumption. Therefore the k — 1st order minimum conditions for / apply: 
for all i = , 7]) G IR2 with f • n > 0, 

i=o \ ^ 



a 



Applying Lemma [3. 1.21 to (|3.1.5p and p.l.6p . with aj = ^ and using 

L < 1, yields D'^V = 0. Note that 6 ^ iff (|3.1.4p is not satisfied. The induction 
step is complete. 

We have shown that D^ip = for all k > 2. Since ip is analytic, it must be 
linear which represents constant density and velocity. □ 

Remark 3.1.4. Lemma 13.1.31 applies trivially to the interior case U: any 
n 7^ will do. 



3.2 Density in the interior 

Proposition 3.2.1. Let x be an analytic solution of \2.2.4^ in an open con- 
nected domain 51. Assume that L < 1 in and that p is positive and not 
constant. Then p does not have maxima in points where Vx 7^ 0, and it does 
not have minima anywhere. 

Remark 3.2.2. Proposition 13 . 2 . II triviallv implies corresponding results for vari- 
ables like p and c that are strictly monotone functions p (except for c in the 
isothermal case where it is constant). 

Proof of Proposition \8.2.1\ The first-order condition for a critical point is 

• V(p) = V(7r-i(-x - ilVxH) = -4V'^VX. 

2 
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If Vx ^ 0, then combined with the PDE (|2.2.5p we obtain D'^^ = 0. Now we 
can apply Lemma l3.1.3l to show that tp is actually a constant-state solution. In 
applying the lemma we choose / = f^ip^Vipya) = f{a) only, using fa^O and 
Vx 7^ so that (|3T4ll is false. 

If Vx = 0, then V{p) = is trivially satisfied. We need to study the second- 
order condition for a minimum, which implies in particular 

< Aip) = (7r-i)'(-Ax - I V^xl^ - Vx • VAx) + (^-^)"| Vx + V^x Vx P 

= pc-2(A?A- IV^VI^) 
The equation (|2.2.5p reduces to A-0 = 0, so 

0<-pc-2|VVr (3.2.1) 

Since p,c> this imphes V^ip = 0. Then V^(/9) = as well. 

We show for fc = 3, 4, 5, . . . by induction that D^{p) = and D^'ijj = as well. 
Induction step (2, . . . , A: — 1 fc > 3): a minimum of p = 7r^^(— x — ||VxP) is 
the same as a maximum of x + ^|Vxp. For j = 0, . . . , k: 

didfix + ^iVxH = didt'x + Vx • Vdidt'x 

^0 

+ jdiVx ■ di-^dt'^X + (fc - .7)52 Vx • didf^^X 
+ terms with components of D^tp, . . . , D'^~^^p as factor 

= (1 - k)didt^iJ 

Here we used that 9iVx = (-1,0), 92 Vx = (0,-1) because V^?/' = 0. The 
induction assumption, D^^p = for j = 2, . . . , fc — 1. eliminates the other terms. 

Since D^{x + ^|VxP) for j l,...,fc — 1, a maximum requires that 
D'^ix + jl^Xp) ^ (i.e. is a negative semidefinite tensor), so > D''{x + 
i|VxP) == (1 — k)D''iJj, so D'^ip > 0. Taking fc — 2 derivatives of the equation 
yields (for Vx = 0) 

did^^-^Ai: ^0 (.7 = 0, . . . , /c - 2). 
Lemma [3.1.11 implies that D^'ip = 0. The induction step is complete. 

Again we have shown that D'^tl' == for all k > 2. Therefore ip, which is analytic, 
must be linear. □ 
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Remark 3.2.3. Note that the proof fails for maxima in V^: in that case < in 
(|3.2.1[) turns into > which does not yield sufficient information. Indeed there 
arc counterexamples. 



3.3 Velocity components in the interior 

Proposition 3.3.1. Let x be an analytic solution of \2.2.4^ in an open con- 
nected domain fi, with p > and L < 1 inil. For any w £ IR^ — {0}, the velocity 
component w ■ V'i/' does not have a maximum or a minimum in fl, unless x "i-s o, 
constant- state solution in Q. 

Proof. Assume that ■01 has a minimum in some point. Then V^i = 0; using 
the equation (|2.2.5p yields D^ip = because L < 1 implies — X2 > 0- Now 
we apply Lemma 13.1.31 to obtain that ip must be a constant-state solution. 

Any other w can be treated by rotating around the origin so that w ■ Vtp becomes 
tjji (see Remark I^XTj) . □ 



3.4 Velocity components on the wall 

Proposition 3.4.1. Consider a point on a straight line I, let r > 0, U := 

Br i^o), r :— I nU and C/+ one of the two connected components of U — I . 
Consider a solution ip of Ili2.2.4\ ) that is analytic in C/+ and satisfies the slip 
condition Xn = onT. Assume that L < 1 in S,o- 

For any w, 

V'0(^o) • w = inf S/ip • w 
is not possible unless ip is a constant- state solution, or unless w is normal to T. 

Proof. Assume there is an extremum point on F. We may assume (by rotation 
and translation) that F is a piece of the horizontal axis, that the extremum point 
is the origin, and that C/+ is contained in the upper halfplane; then w = (w^, w^) 
with ^ (not normal). A tangential derivative of the boundary condition 
tp2 = implies 7/112 = on F. A minimum requires 

= (wVV')l — W^tpii + w'^tpi2 = W^^ll "011 = 0. 

The equation (|2.2.5p yields that D'^tl' == 0. Now the result is delivered by Lemma 
[3T31 □ 



29 



3.5 Velocity components at shocks 

Proposition 3.5.1. Consider disjoint open connected domains fl and $7'* and 
a simple analytic curve S d Q (1 Q (excluding the endpoints). Consider a 
constant- state (linear) potential ip'^ of i2.3.1]) in fi''. Let x be an analytic solu- 
tion of {2.2-4^ in flu S. Let the shock relations f2.3.1\) and Ii2. 3. 2\) he satisfied 
on S and assume the shock is admissible. Assume that x satisfies p > and 
L<1 inflU S. 

Let w ^ 0. Assume that w ■ Vi/' has a local maximum (with respect to flU S) in 
^ G S". Then either S is straight and ip is constant- state in D,, or 

w- ({l-c-\l)t + xti^ + c-\n)n) =0. (3.5.1) 

Xn 

and 

sgnK = sgnii;" ^ (3.5.2) 
where k is the curvature of S in ^ (k > for locally convex). 

Proof. We use a dot to indicate relations that hold only in the hypothetical 
extrcnium point. By L < 1, fl must be downstream and fl^ is upstream. 
Without loss of generality, rotate around until n = (0,-1). In this setting 
di = dt and —82 = dn, and xi ^ Xt — Xt ^ x'i by (|2.3.4p . Use horizontal 
translation (Remark 12.2. ip so that ipi = and therefore ipi = 0; this adds a 
constant vector to velocities, while leaving p and L unchanged, so no generality 
is lost. Let the shock be parametrized by ^ i-^ (^, s(^)) locally; then si = 0. 

Take du of (^X^ : 

^11 + (V'2 - ^2 )sii = 0. (3.5.3) 

Take dt of (l2X2l) : 

= dt {{tt-H-x - ^IVxHvx - /vx'^) • n) 

= n ■ dt{7T-\-x ^iVx^Vx - p'^^x") + (t-^(-X - ^iVx^Vx - /Vx'^) • 
• 9i(^-i(-x - ^IVxHVx - /Vx'') + [tt-H-X - ^iVxHVx - p'^Vx 

= pc"2(xi + XiXii + X2Xi2)X2 - PX12 + p''Xi2 + (PXi - p''Xi)sii 

= p( - (1 - c-\l)^u + c-\iX2^ii + Xi(l - (3.5.4) 
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Combining these results with the equation and {w ■ V'0)t = we get the system 

,1 -,,,2 



w w 
1 

- xl -2xiX2 

,-2,, 1 „-2,,2 







? 9 
C -X2 



' X2 








V'li 


-x^ 




V'12 







V'22 


-i)xi 




Sll 



0. 



(3.5.5) 



_-c ^XiX2 1 - c ^X2 
Determinant of the system matrix: 

det = -(c2 - x^)(^2Xi _ _ c-^x^) + u;2c-2xiX2)) 



X2 

:(c'_^^x^)(x2_^x^'W 

>0 >0 



1 - c-^xl 

Xl(l/X2 +C-'X2) 



(3.5.6) 



The determinant is zero iff the final scalar product is zero. The latter condition 
can be written (|3.5.1[) . if we return to original coordinates. 

If the determinant is nonzero, then Z)^ V ^i^d sn = is the only solution. If 
the determinant is zero, but sn = 0, still ~ 0. In either case we can invoke 
Lemma [3. 1.31 to get that is a constant-state solution; then (|2.3.3p shows that 
the shock is straight. 

Now assume the determinant is zero and sn ^ 0. By row 2 of the system this 
implies V'li 7^ 0. Then by row 1 and it; 7^ 0, necessarily vP' 7^ 0. 

By solving rows 1,2,3 of the system for V^V as a function of sn, then substi- 
tuting the result into (w ■ VV')2, we obtain 



(w ■ VV')s 



(X2 - X2)((c^ ~ xi)wi + 2xiX2W^w^ + (<? - Xi)wi) 

W2(c2 _ 

{X2-xfiKw^nc^I-^x')w^ sn 

C2 - xl W2 



Sll 



All factors in the coefficient of Sii/w2 are positive; sn, W2 7^ is already known, 

so {w ■ V?/')2 7^ 0. A maximum requires {w ■ V?/')2 > 0, so sguK = — sgnsn = 
— sgnit;^ = sgntii" 7^ 0. □ 



3.6 Pseudo-Mach number at shocks 

Proposition 3.6.1. Consider the setting in the first paragraph of the statement 
of Proposition \3.5.1\ 

^Wc do not need analyticity here. 
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Let > be such that 



— e [Slp. 1 - Slp\ 



s 



(3.6.1) 



Let b e C^{flU S). There is a Si^g > (depending continuously and only on 
6lp,J, c^, p'V "with the following property: 

+ b cannot attain a local ( with respect to S) maximum in a point on S 
where e [1 — 6ls, 1) and |V6| < Sls- 



Proof. We use the same notation and simplifications as explained at the start 
of the proof of Proposition 13.5.11 



From pI < 1 — dLp < 1 we obtain 



C2 ' X2 



(3.6.2) 



for some constant Cg = Cs{5Lp,"i) > 0. 
diiL^ + b) 



(3.6.3) 



-) + bi 



^ci^ + (l-7)(x+i|VxP)- 
= c-2((2 + (7 - l)L2)(xiXii + X2X12) + (7 - l)i'xi) 
= (^(2 + (7 - 1)l2)(xiV'ii + X2V'i2) - 2xi) + bi 
and analogously 

d2iL^ + fe) = ((2 + (7 - 1)L2)(xi^12 + X2V'22) - 2X2) + &2 



bi 



(3.6.4) 



Combining 5i(L^) = with the equation (|2.2.5p . and p.5.3p and (|3.5.4p we 
have a linear system 



Xi 


X2 





1 








9 9 

- XI 


-2xiX2 


9 9 

c-xi 


XlX2 


1-4 








X2 - X2 



(fl-i)xi^ 



i'll 
i'i2 

i'22 
Sll 



2+(7-l)L^ 







First consider V6 = and L ~ 1, i.e. xi = sgnxi • ■\/c^ ^ X2 • Then the inverse 
of the system matrix has entries polynomial in c, X2, X2 1 ^g^iXi divided by a 
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common denominator 



This denominator is bounded below away from zero by c' 5d where 5d depends 
continuously and only on S^p and 7. The rest of the inverse matrix is bounded 

Lp- 



by some constant depending only on p'', c'\ 7, 5 



Solving for D'^^jj and substituting the result yields 

X2 + X2 

SO clearly a maximum of + 6 is not possible. 

For V6 ^ and L < 1 wc use that the inverse has been bounded away from 
0, so that small perturbations are possible. Thus, there is a Sls, depending 
only on Slp, 7, p'^, ■, so that no maximum is possible if |V6| < 5ls and 
> 1 - Sls- □ 



Remark 3.6.2. We could assume L < L < 1 which by itself would imply 
p/p^ < C{L) < 1. However, this is not sufficient as 6ls would depend on 
L, so Proposition 13 . 6 . II would be void. It is necessary to obtain uniform shock 
strength bounds separately; only then can L be controlled. 



3.7 Density at shocks 

Proposition 3.7.1. Consider the setting in the first paragraph of Proposition 

If p has a local extremum with respect to flU S in ^0, then one of the following 
alternatives must hold: 



1. X a constant- state solution in QU S , and S is straight. 

2. The shock is pseudo-normal in ^o- In a local minimum, S has curvature 
K > (as before, k > for fl locally strictly convex in that point). 

In the latter case, assume stronger that p has a global minimum with respect to 
$7 in S,o- Assume for technical convenience that the shock tangents differ by no 
more than an angle < 7r/2 from the one in ^q. Let S* be the tangent to S in 
i^Oi s{Co))- Then S* does not meet S anywhere else. 

Proof. We use the same notation and simplifications as explained at the start 
of the proof of Proposition 13.5. II 
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A p extremum requires 



= dtip) = 9i(7r-i(-x - ^iVxH) - -^{Xi + XiXii + X2X12) 



Xl-iAll + X2^'12 = 



We combine this with the now-famihar equations (|2.2.5p . p.5.3p and (|3.5.4p . 
The resulting hnear system is 



= 0. 



Xi 


X2 













C^-Xl 


-2xiX2 


- X2 


-c ^XlX2 


i-c-^xi 



















-x'i 




V'12 







V'22 


-i)xi 







The determinant is 



-(c'"X2)(X2-X^)Vx 



i-c-^xi 



XI 
X2. 



c "X1X2 



-(c^;^_xg(x^^(l + ^)-Xi 
>o >o V ' 

>0 



It is nonzero if and only if xi 7^ 0. In that case, V^V = a-^d sn = 0. Now 
Lemma [3.1.31 yields that the shock is straight and the solution constant-state. 
(The Lemma applies because p is a strictly decreasing function of x + ^I^Xp 
alone, and p.l.4p is satisfied because Xn 7^ 0, hence Vx 7^ 0, at any shock). 



If Xi == but sii = 0, then the equations imply that V^?/; 
still applies. This concludes the proof of the first part. 



0, so the Lemma 



For part two we may assume that all of S is parametrized by s, because the 
shock tangents cannot become vertical, by assumption. Consider the remaining 

• * • • • 

case Xi = ^'iid sn ^0. Xi = ^'iid si = imply ^ = 0. Here it is sufRcicnt to 

argue without the interior. So we may exploit that at the shock downstream, 
Pd is an increasing function of z" (by (|2.4.10p and (|2.4.13p . for p„ held fixed). 
Geometrically, 2" in a point ^ on the shock is the distance of Vu to the shock 
tangent through ^. A pseudo-normal point is actually the closest point on the 
tangent to Vu- Obviously the tangents through nearby points are closer to u„ 

if sii > (sec Figure ini), which contradicts the assumption that p is a local (in 
fact global) minimum. Therefore sn < 0. sn = has already been excluded, 
so fill < 0. 

Let ^ ^ s*(0 parametrize S. Assume that S* meets S somewhere else, e.g. that 
s(^i) = s*(^i) for some ^1 < ^0 (see Figure ITO)) . By the mean value theorem 
there must be a ^2 € (Ci,$o) so that 51(^2) = (s(6) - s(^o))/(Ci - Co) = si($o)- 
Clearly we may choose ^2 maximal with this property, so that si(^) ^ si(^o) 
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sii > 



Shock 




Shock lower, 
parallel tangent 

higher pa 
Contradiction! 



Global pd minimum 
sii < 




Figure 9: Tangent distance ar- 
gument 



Conclusion: impossible 

Left corner must be below min point tangent 

Figure 10: Shock below global density mini- 
mum tangent 



?/)-normal point 




Figure 11: Shock strength: distance of tangent from upstream velocity 
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for C € (6,6)- But fill < 0, so si(^) > si(^o) for C S (6,^o)- This implies 
5(^2) < s*(Co): the shock tangent through (^2,5(^2)) is parallel to the one in 
^0, but lower. That means the shock strength is smaller, so the downstream 
density is lower (by ()2.4.20|) ) — contradiction, because we assumed that p has 
a global minimum in ^o- 

Therefore s(^) > s(^o) for all ^ < (on S). For ^ > ^0 the arguments are 
symmetric. □ 



4 Construction of the flow 



4.1 Problems 

Our solution, as observed in numerics (see Figure[4|), has the structure in Figure 
[T2I The upstream region, labeled is a constant-state hyperbolic region. 
The shock has three parts: a straight shock emanating from the tip, with a 
constant-state hyperbolic region labeled "L" below; a curved shock S with a 
nontrivial elliptic region below; and a straight shock parallel to the wall, with 
another constant-state hyperbolic region ("i?") below. The L and R regions are 
separated from the elliptic region by parabolic arcs Pl resp. Pr with radius cl 
resp. cr, centered in vl resp. vr. 

Several difficulties complicate the problem: the equation is nonlinear (quasilin- 
ear, divergence form), with coefficients depending on x o-^d Vx- The boundary 
conditions (except on the wall) are fully nonlinear, i.e. they are nonlinear in Vx 
as well, which makes compactness hard to obtain. Moreover, the boundary con- 
ditions linearize to oblique derivative conditions where the x and Xn coefficients 
have opposite signs — the most complicated case. The equation and boundary 
conditions have singularities: for example we have to avoid the vacuum {p ~ 0). 
The shock is a free boundary, forming twcj^ corners with the arcs. 

Most significantly, the equation is mixed-type: if x and Vx are not sufficiently 
controlled, points in a supposedly elliptic region could be parabolic or hyper- 
bolic. While the elliptic region is uniformly elliptic at the shock, it is degenerate 
elliptic at the parabolic arcs P^ and Pr. Moreover it appears from numerics 
that Vx is normal on the parabolic arcs, meaning they are characteristic (in 
the Cauchy-Kovalevskaya sense), so loss of regularity has to be expected. The 
linearization of the degenerate problem is useless in this case. E.g. although we 
expect finite gradients {— velocities) in the nonlinear problem (and observe them 
in numerics), it can be checked for simple examples that the linear equation has 
solutions with gradient that is infinite at parabolic arcs. 

For many of these obstacles there are theoretical tools in the literature; many 

^The wall-arc corners can be removed by reflection across the wall. 
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have been addressed in other contexts. However, the nonhnear characteristic 
degeneracy seems unprecedented. There is little theory on degenerate elliptic 
equations; most of it can be found in [UR73j . A large part of the theory was 
motivated by the linear Tricomi equation which arises from steady potential 
flow via the hodograph transform (see |Eva98[ Section 4. 4. 3. a]). The hodo- 
graph transform applies to quasilinear equations whose non-divergence form 
coefficients depend only on the gradient of the solution; the original equation 
is converted into a linear equation. Unfortunately the coefficients of selfsimi- 
lar potential flow also depend on the solution itself. It is not clear whether a 
modified hodograph transform can be developed for selfsimilar potential flow. 

Nonlinear degenerate elliptic equations, other than steady potential flow, have 
not been explored much, apart from steady potential flow which can often be 
reduced to a linear equation. The combination of nonlincarity and degeneracy 
is particularly difficult: the solution can be characteristic degenerate, nonchar- 
acteristic degenerate, elliptic or hyperbolic in each point of the arcs; each case 
is qualitatively very different. Precise bounds on the solution have to be estab- 
lished before we even know which case occurs in which point. 

4.2 Approach 

In this paper we construct a weak solution, so only x continuity and conservation 
are shown across the parabolic arcs. Most of the effort is concerned with the 
elliptic region. 

Symmetries Due to Remark 12.2. 1[ we can change our coordinate system in 
several ways. 

Starting in original coordinates (see Figure 5]) , we translate ^ — vji and 
V ^ V ~ vr, then rotate so that the right shock is horizontal (and in the 
upper half-plane); see Figure [T2l Now vr = 0, vl = {vf^,0) with vf^ < 0, 
and vi = {0,Vj) with < 0. Moreover, ip is constant in the R region, and Pr 
is centered in the origin, which will be rather convenient. We use this as the 
standard picture, as it is the most convenient for discussing the elliptic region, 
which consumes most of our efforts. 

For another choice, which we call "L picture" , we start in original coordinates, 
translate ^ <— — and v ^ v — v^, then rotate so that the left shock 
is horizontal (and in the upper half-plane), then reflect everything across the 
vertical coordinate axis; see Figure [T^l In this picture vl = whereas vr = 
{v^, 0) with vfj < 0, and vj = (0, ) with < (which need not have the 
same value as before). Now ip is constant in the L region, and Pl is centered in 
the origin. 
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Figure 12: Left: standard coordinates; right: L picture 



Extremum principles As for most nonlinear elliptic problems, maximum 
principles are the key technique. In Section [3] we have developed many that 
apply to selfsimilar potential flow in general. Their interior versions are compa- 
rable to the classical strong maximum principle |GT83| Theorem 3.5] or to max- 
imum principles for gradients of special quasilinear equations [GT83[ Section 
15.1]. In addition we have several extremum principles at the shock (Propo- 
sitions I3.7.1i 13.5.11 and 13.6. 1[) , ruling out local (with respect to the domain) 
extrema of certain variables at the shock. 

The general proof technique for extremum principles is to combine the equation 
with the first and second order conditions for an extremum to obtain a contra- 
diction. At the shock we also include the boundary conditions. In many cases 
it is necessary to include first and higher derivatives of equation and bound- 
ary conditions as well as higher-order extremum conditions. In some borderline 
cases (notably density and velocity) it is necessary to consider all derivatives 
and to assume that the solution is analytic (see the proofs of Propositions 13. 2. II 
[HTmiOTTI and lXCT . 



Regular izat ion and ellipticity As we have mentioned, our elliptic region 
is nonlinear and characteristic degenerate; there are very few results for such 
problems. The type of our PDE is governed by the pseudo-Mach number L. 
We avoid the degeneracy by regularization: we consider slightly different elliptic 
regions, with L = VI ~ e on the parabolic arcs rather than L ~ 1 (see Figure 
I13p . We need L uniformly bounded above away from 1, so that our regularized 
problems are uniformly elliptic and standard theory can be applied. 

A basic problem of our constructive approach is that it is infeasible if the flow 
pattern has a complicated structure. For example if the elliptic region could, 
upon perturbation, develop a (potentially infinite) number of parabolic and 
hyperbolic bubbles in its interior, only abstract functional-analytic methods 
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Figure 13: The solution in the eUiptic region and its dchmiting shock S are un- 
known. The paraboUc arcs (sohd) are modified to quasi-parabohc arcs (dashed) 
with boundary condition L = \/l — e which is slightly elliptic. 

have a chance to succeed. Fortunately we have shown in |EL05a| that the 
pseudo-Mach number L cannot have maxima in the interior of an elliptic region 
or at a straight wall. On the arcs we have imposed the value of L. It remains 
to control L at the shock. This result (Proposition [21111]) is less general: we can 
rule out L maxima close to 1 at the shock if the shock has uniform strength. 
This condition needs to be verified with separate methods (see below). 

All combined we can be sure that our elliptic region stays uniformly elliptic 
under perturbation. We obtain, after many additional steps described below, 
one solution x^^-* for each sufficiently small e > 0. To obtain a solution of the 
original degenerate problem it is necessary to obtain estimates that are uniform 
in e i 0. 

As discussed in jEL05a| , the maximum principle for L can be strengthened to a 
maximum principle for L + b where & is a small smooth positive function which 
is zero on the arcs. As a consequence, L is uniformly m e J, bounded above 
away from 1, in subsets of the elliptic region that have positive distance from 
the arcs. We obtain uniform regularity in all C*''" norms (in fact analyticity) in 
each such subset. Hence we have compactness in each subset, so we can find a 
converging subsequence. By a diagonalization argument we have a subsequence 
that converges everywhere away from the arcs, in any norm (see Section [4. 16p . 

Another of the many benefits of the ellipticity principle: L < 1 ^ e < 1 implies 
a uniform bound on the gradient |Vx| (see (|4.5.5|) etc), which is a basic ingre- 
dient of all higher regularity estimates for nonlinear elliptic equations (see the 
discussion in [GT83[ Section 11.3]). Moreover it shows that x^'' is uniformly 
Lipschitz, so x'-"^ is Lipschitz as well, including the arcs. This is not quite suf- 
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ficient for the boundary condition i = 1 to be well-defined; but we only seek a 
weak solution. 



Iteration For any nonlinear elliptic problem, there are several ways of turning 
a priori estimates into an existence proof (cf. [GT83| Chapter 11, Section 17.2]). 
The method of continuity is straightforward: verify that the linearizations of the 
problem are isomorphisms between suitable spaces, apply the inverse function 
theorem to obtain small perturbations, then use the a priori estimates to show 
that we can repeat small perturbations indefinitely. 

A major drawback is that the linearizations are entirely determined by the non- 
linear problem. In comparison, fixed-point iteration methods are more flexible: 
there is a large variety of maps whose fixed points solve the nonlinear problem. 
In our case this flexibility is vital: the Frechet derivative of = 1 — e is an 
oblique derivative condition where Sx and Sxn i^X being the variation of x) 
have opposite sign. In this case, maximum principles fail, but they are needed 
to verify uniqueness (the other option, energy methods, seems unpractical for 
our complicated problem). 

Another (non-fatal) drawback is that linearizations of a free-boundary problem 
involve a complicated coordinate transform. Iteration methods can alternate 
solving a fix ed-boundary problem with adjusting the free boundary (see e.g. 



ICanicKLOO) and ICFOSQ . which is easier. 



Fixed-point methods can be subdivided into applications of the Schauder fixed 
point theorem (or its generalizations), and of Leray-Schauder degree theory (see 
|Dei85j . [Zei86] or (SmoSl)- 

The Schauder fixed point theorem requires showing that the iteration is maps 
a closed ball (or homeomorphic image thereof) into itself. This corresponds 
to certain estimates for all elements of the ball, most of which are not fixed 
points (only approximately perhaps). These new estimates are particular to the 
chosen iteration, which is unattractive and perhaps difhculll!!. Moreover, every 
time details of the iteration are changed, all the a priori estimates have to be 
checked and changed; this takes an excessive amount of time for complicated 
problems like the present one. 

In addition, while many function sets in nonlinear elliptic problems are defined 
by a single constraint ||u||,72.q < M < oo, so that the set is obviously a closed 
ball, our function set J- is defined by more than twenty different constraints 
(see Definition I4.4.3p . For many choices T has nontrivial topology; it is rather 
difficult to check whether an intersection of many sets is homeomorphic to a 
ball (unless convexity can be used). 



^However, ICF03I have successfully applied the Schauder fixed point theorem to solve a 
steady potential flow problem. 
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Ultimately the Schauder fixed point theorem is a means of showing that a par- 
ticular iteration has nonzero Leray-Schauder degree. However, there are other 
ways of computing a degree which turn out to be simpler in our case (see Section 



The iteration has two steps. The first step solves a fixed-boundary elliptic 
problem for a function hi this step we use a modification of the second 
shock condition (I2.3.2[) : the data of this problem depends on x (argument of the 
iteration map). In the second step the shock is adjusted (and x mapped to a 
new domain) to satisfy the first shock condition (|2.3.ip . 

As mentioned, we have some flexibility in choosing the first-step elliptic problem. 
An obvious constraint is that fixed points of the iteration, x = X^ must solve 
the original problem. This is satisfied by replacing every occurence of x iii the 
original equation and boundary conditions either by x (old iterate) or by x (new 
iterate after step 1). 

For the arc boundary condition (|4.7.ip we choose 



for the iteration, where x is the old and x the new solution. Linearization with 
respect to x no longer has a 0th order term! Hence the opposite-sign problem 
does not apply, and we can apply the Hopf lemma to show that the linearizations 
of the nonlinear elliptic problems arising in the iteration are isomorphisms. We 
use the same technique for the interior equation, but keep x in ttic p function 
in the shock condition (|2.3.2p : we need at least one occurence of x to be able to 
use maximum principles to achieve uniqueness in various contexts. Obviously if 
X does not appear anywhere in the elliptic problem defining the iteration, then 
the next iterate is not uniquely determined: any constant can be added. 

To apply Leray-Schauder degree theory, we argue that we are solving a contin- 
uum of problems, the simplest one being the "unperturbed" problem (see Figure 
I^O)) with 7=1. This problem is simple enough to verify that its Leray-Schauder 
degree is ±1, in particular nonzero (see Section l4.14p . All other problems have 
the same degree, due to the combined a priori estimates (no fixed points on the 
boundary of the iteration domain) and continuity. 

The opposite-sign difficulty resurfaces in two aspects: first, to prove degree 
7^ we have to show uniqueness of the unperturbed problem (see Proposition 
I4.14.ip . For 7 > 1 this seems difficult, but for 7=1 the x coefficient is zero (see 
above). Moreover, it is necessary to show that the linearization of the iteration 
does not have eigenvalue 1. Again 7 > 1 has unpleasant boundary conditions, 
but 7 = 1 is amenable. Thus solving the isothermal and isentropic problems in 
a single paper is necessity rather than choice. 

Unfortunately it is necessary to show that the iteration is compact, which re- 
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quires somewhat stronger regularity estimates than the method of continuity. 
In fact |GT83[ page 482] state: "Even for the quasihnear [equation] case [with 
fully nonlinear boundary conditions], the fixed point methods [...] are not ap- 
propriate, since it is not in general possible to construct a compact operator 
[...]." 

Recent results, many of which are due to Gary Lieberman, have alleviated this 
problem, at least for our case. Some unpleasant choices are necessary, however: 
by replacing the occurences of Vx by a clever balance of Vx and Vx we could 
obtain a linear elliptic problem as part the iteration. But if Vx occurs in any 
of the boundary conditions, then the new solution x cannot be expected to be 
smoother than x at the boundary. While we may use Xi Vx in the interior 
coefficients, we must use Vx for the boundary conditions, which stay nonlinear. 

Now our general approach has a hierarchy of three elliptic problems: the original 
nonlinear problem, a modified nonlinear problem as part of the iteration, and its 
linearization. It appears that we have to start all over showing a priori estimates 
for the modified nonlinear problem. A triclH avoids this: if the linearizations 
evaluated at x = X are isomorphisms, then the nonlinear problems are local 
isomorphisms. Thus we restrict the set of x so that the next iterate x = ^(x) is 
close to x- Existence and uniqueness of x as well as its continuous dependence 
on X then follow by linearization around x; large perturbations are not necessary. 

Boundary of the function set Both method of continuity and degree theory 
have a common feature: the bulk of the effort goes in showing that the problem 
does not have solutions on the boundary of the function set. For function sets 
defined by inequality constraints, we have to show that for a solution of the 
original problem (a fixed point of the iteration), each inequality is in fact a 
strict (<, >) inequality. In this step we are allowed to use the nonstrict (<, >) 
forms: wc are moving from one well-behaved (smooth, elliptic, no vacuum, 
...) solution to another, without having to show any results for arbitrary weak 
solutions with no prior information. In particular wc do not have to verify the 
strict inequalities in any particular order. 

Parabolic arcs and corners The core of the solution — and its most difhcult 
part — is the treatment of the parabolic arcs: to obtain a weak solution in the 
limit, we have to verify the shock conditions are satisfied across the arcs. To this 
end we show that x^'^ and Vx^*^^ are almost (up to 0(e-^/^) jumps) continuous 
across the arcs; in the e J, limit continuity and mass conservation for x^"' are 
implied. 

* Another, less elegant, trick is to use the a priori estimates of the original problem to cut off 
the coefficients of the second problem so that they grow linearly and avoid all singularities like 
vacuum or parabolicity. The cutoff terms and their derivatives would make the linearization 
unnecessarily complicated. 
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X, Vx continuity corresponds to two boundary conditions, 

lim X = lim Xi 1™ Xn = lim Xn 

hyperbolic elliptic hyperbolic elliptic 

(continuity for xt is implied by continuity for x)- However, with two boundary 
conditions the problem is overdetermined, at least in the regularized case. Hence 
we can impose only one boundary condition. A trick is needed to verify the 
second condition — approximately for the regularized solutions; exactly (albeit 
weakly) in the degenerate limit: 

In each point on the arcs, there are three components of V^Xi but only two 
relations constraining them — the equation (j2.2.4[) and the tangential derivative 
of the boundary condition, = 1 — e (see (|4.7.2|) ). Our key observation is that 
< 1 — e inside the elliptic region (as we have already shown previously), 
whereas = 1 — e on the arcs: L has a global maximum in each arc point. 
Therefore (i^)n > on the arcs, which provides an additional inequality. Thus 
we obtain "2^" relations; solving them yields one explicit inequality for each 
second derivative, with right-hand side depending only on x f^nd Vx- Only the 
twice tangential derivative matters (see (|4.7.6p ). 

In the shock-arc corners x is (Proposition 14.5.^ using Section [5TT|) . so we 
may combine the boundary conditions (one on P, two on the shock, but with 
unknown shock tangent) to represent x and Vx as functions of the corner loca- 
tion. If we assume that the corners have precisely the "expected" location, then 
we know the corner values of x and Vx; in particular xt = 0, where dt is the 
counterclockwise tangential derivative along P, and c = (upper left corner) 
resp. c = Ch (upper right corner). However, our analysis is greatly complicated 
by having free corners. In the other endpoint of each arc the wall boundary 
condition Xri = fixes Yt = 0. Then the ordinary differential inequalities men- 
tioned above yield tighv^ bounds for x and Vx on the arcs: they have the desired 
values, up to 0(e). This fortunate circumstance is the key to our solution. 

We first discuss isothermal flow which is comparatively simple. In this case the 
inequality for xt on has the form 

Xtt >/(Xt) = -e + 0(x?), (4.2.1) 

(see ()4.8.7p ). with O term constant independent of e. Xt ~ in the wall-arc 
corner, so integration to the arc-shock corner yields 

Xt > -0{e) (4.2.2) 

along the entire arc. 

In Section [4.91 we compute the values of x, Vx and dependent quantities in the 
right shock-arc corner as the corner location varies along the arc. (|4.9.10p shows 

^It is worth noting that if, for whatever reason, we had > a for some bounded 

(smooth) function a, the trick would still work. The tightness does not stem from a = (in 
fact is far from optimal), but rather from the characteristic degeneracy. 
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that Xt decreases uniformly as the corner moves down from its expected location. 
This contradicts Xt > —0{e) — the second- derivative inequalities already imply 
(see Proposition l4.10.T|) a lower bound for each corner, to within 0(e^/^) of their 
expected location. 

However, an entirely different argument (see Proposition I4.10.6p is needed to 
bound the right shock-arc corner above: tp niust attain a global minimum in 
the domain. It can never attain a local minimum in the interior, on the wall, 
or on the left arc. For isothermal flow we have -i/in > on the right arc, so the 
Hopf lemma excludes a minimum there as well. Hence there has to be a global 
minimum on the shock or in the right shock-arc corner. Essential observation: 
if that corner is above its expected location, then -02 > in it, as well as in any 
hypothetical global minimum on the shock, so wc obtain a contradiction. 

Having obtained upper and lower bounds on the corner, we know that xt = 
0(e^/^) in it. Integrating (|4.2.1I1 again, but in opposite direction, we obtain 

Xt < -0{e'/') (4.2.3) 

Combined with (|4.2.2[) we have xt = 0{e^^^) on the arcs. By integration we 
also control x; the boundary condition |VxP = yields control over Xn- 

The arguments for the left arc are analogous, with a few modifications. 

For isentropic flow the arguments are similar, but much more complicated, 
because the sound speed c can vary. Instead of considering a single-variable 
ordinary differential inequality we have to combine it with an ODE for c (see 
(|4.7.7p ). To make the system more tractable it is linearized ( (|4.7.10p . ()4.7.1ip ) 
and restated in polar variables ( (|4.7.12|) . ()4.7.13|) ). A delicate analysis (see the 
proofs of Propositions 14.8.11 and I4.10.ip again establishes a lower bound for the 
right shock-arc corner, as well as a lower bound for xt- For an upper bound 
on the corner we adapt the isothermal argument, considering minima of Tp + 
instead of ip, for some small a > 0. This is necessary because for isentropic flow 
i'n ^ —0{e^^^) rather than > on the right arc. More delicate analysis (see 
the proof of Proposition 14.10.61 and its references) shows that a can be chosen 
so large that {tp + a^)„ > 0, but so small that still {tp + a^)^ > in the right 
shock-arc corner and in any hypothetical minimum point on the shock. Having 
bounded the corner above, an upper bound for xt on the arc is obtained as well. 

An entirely different approach to parabolic arcs can be found in [CFj . 

Shock control It is clear that we have to expect degeneracy at the parabolic 
arcs. However, there can be degeneracy on the elliptic side of a hyperbolic- 
elliptic shock, too, which we have to rule out. Moreover we need some control 
over the shock location and shape. 

To prove that the shock has uniform strength, we show that the density on the 
elliptic side is uniformly bounded below away from pi . Proposition 13.2.11 rules 
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out density minima in the interior of the elhptic region or at the wall (Remark 
I2.2.ip . The analysis described above controls the density on the arc up to 
0(e^/^) (sec (|4.8.1|) '). Finally, Proposition 13 . 7 . Il shows that density cannot have 
local minima at the shock except in a pseudo-normal point. The shock curvature 
must be positive (elliptic region locally convex) in such a point. Moreover if we 
have a global minimum in this point, then the rest of the shock, including the 
shock-arc corners, must be below the shock tangent in that point: otherwise 
there would be another shock tangent parallel to this one, but lower, meaning 
lower density, which is a contradiction (see Proposition l4. 1 1 .11 Figure[Tn]for the 
detailed argument). 

Hence we can have density minima at the shock, but only in pseudo-normal 
points and above the line connecting the shock-arc corners. In such a point, the 
density is > pi if and only if the connecting line docs not meet the circle with 
center vi and radius c/. (Otherwise the shock vanishes or is entropy-violating.) 
This is precisely (jl.3.ip . Unfortunately the cases where (jl.B.ip is violated are 
not covered in the present paper. 

Having established uniform shock strength, Proposition l3. 6.11 shows that L can- 
not have local maxima (with respect to the domain) close to 1 at the shock. This 
is a necessary ingredient for L control in the elliptic region; the ellipticity must 
be uniform at the shock (away from the corners). Hence the shock is analytic, 
except perhaps in the corners. 

To control the shock tangents, p and L arguments arc not sufficient. It is 
necessary to control some components of the velocity vector. More precisely we 
control the horizontal (in standard coordinates) velocity u^, as well as horizontal 
velocity in the L picture which is + av^ in the standard picture, for some 
a > 0. We show that v'^ cannot have maxima in the interior fProposition l3.3.l|) . 
at the wall fProposition l3.4.l|) . or at the shock (Proposition 13. 5. l|l . In the latter 
case, like for density, there are exceptional cases (see (|3.5.ip ) where the shock 
curvature (see p.5.2|) ) is needed to rule out maxima. On the arcs we control 
velocity (like everything else) up to 0(e^/^) (see (j4.8.2p ). 

All combined we have that must be between vf^ < and vfj = 0, up to 
0(e^/^), and an analogous result in L coordinates (see (|4.4.24p . (|4.4.26|1 ). This 
yields a slew of additional information (Proposition l4.12TT|) : bounds on the shock 
normals (|4.4.27p . the distance of shock to wall ()4.4.4|) . the shock-arc corner 
angles (j4.4.6p . and the vertical velocity (|4.4.25|) . 

Note that the vertical velocity can have minima at the shock, as can be observed 
in numerics (Figure [S] bottom) . Other velocity directions can have extrema as 
well; we would be able to cover all cases of the theorem, even those violating 
(|1.3.ip , if we had perfect velocity control. 

As a convenient sideeffect we have ruled out vacuum or negative densities, which 
is another type of singularity affecting self-similar potential flow. 
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4.3 Parameter set 



We consider 7 G [l,oo), pj E (0,oo), cj e (0,oo). In addition we use Mj G 
(—00,0) which defines vj = {0,Mfci) (in standard coordinates; see Figure [T^ . 
Mf is the upstream Mach number normal to the downstream wall, pj and vj 
define a potential for the / region: 

Given upstream data Vu = vi and Pu = Pi, Proposition l2 .6 . 21 defines a horizontal 
shock {R shock) with vn := Vd = 0. Let ?7|j be its height; we choose •q*^^ as a 
parameter that determines u/, rather than vice versa, so that we may regard ry^ 
defined as independent of 7 (of course w/, depend on 7 now). 

Let e G (0,e] (for some suitable small e > which will be fixed later). Of the 
two Ld = Vl — e points for this shock, as defined by Proposition 12 . 6 . ll let 
be the one farther from the origin (see Figure [72]) . Set fifi = (0,-1). 

Starting in the R shock. Proposition 12 . 6 . 2] vields a family of shocks with = 0. 
Each has two Ld = \^T~~e points; let ^2 be the one closer to the origin. We 
focus on choices rjl G (0, ry|j]. We call this shock the L shock. Let ul he the 
corresponding downstream normal. ^2 and will be called the expected corner 
locations (although they are most likely not the true locations, except for e = 0). 

Let cq, Pcj vc be the downstream data of the C shock (C G {L, R}). Note that 
vl = 0, but vl <0 for ?y2 < Vr- 

Define ^br := {cr,0) and £^bl = (w£ - cl,0). Let A = {£.bl,^br) x {0}; we 
will call A the wall (it is only the elliptic portion of the wall). Let be a 
circular arc centered in vn passing counterclockwise from ^br to let P£ be 
the arc centered in vl passing counterclockwise from ^2 to ^bl', both excluding 
the endpoints. Pc has radius ^/T~~e ■ cc (for C ^ L, R). 

There is a Snt > so that, for any unit vector ts from nj^ to nj^ (counterclock- 
wise) and any unit tangent tp of P£ or P^, 

\ts X tp\ > Snt > 0- (4.3.1) 

We choose extended arcs Pl,r that overshoot ^2 ^ by an angle Sp > 0, which 
we choose continuous in "f,r]1, so that 

\ts X TpI > 6nt/2 (4.3.2) 

for the same ts, but unit tangents tp of the extended arcs Pl,r- Le. the possible 
shock tangents (restricted in (|4.4.27p below) and arc tangents are uniformly not 
collinear. 
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P£ fli Pl,fi, and later Pl,r, are called quasi-parabolic arcs (or parabolic arcs, 
by abuse of terminology, or short arcs). (Of course these arcs are circular; 
"parabolic" refers to the expected type of the PDE at these arcs.) 

Wc use Pi/j^ and P£ 'j^ to identify the arcs for a particular choice of e. 

The Dcfinitions l4.3.1[|4.4.2l and l4.4.3l use many constants and other objects that 
will be fixed later on. In all of these cases, an upper (or lower) bound for each 
constant is found. Whenever we say "for sufficiently small constants" (etc.), we 
mean that bounds for them are adjusted. To rule out circularity, it is necessary 
to specify which bounds may depend on the values of which other bounds. In 
the following list, bounds on a constant may only depend on bounds of other 
constants before them. 

5o,Cd,e,Cc,rj,a,(3. (4.3.3) 

The constants Cc,ri,a, (3 may depend on e itself, not just on an upper bound, 
r/ may also depend on ip. The reader may convince himself that the following 
discussion does respect this order. 

The parameters 7, jy^ used in Leray-Schauder degree arguments will be re- 
stricted to compact sets below so that any constant that can be chosen continu- 
ous in them might as well be taken independent of them. Dependence on other 
parameters like pi will not be pointed out explicitly. 

Constants S-/ as well as a, /3— 1, r/, e are meant to be small and positive, constants 
C? are meant to be large and finite. 

Definition 4.3.1. For the purposes of degree theory we define a restricted 
parameter set (sec Figure IM)) 



with 

where 



A := {-^ = (7,'/l) : 7 e [1,7], vil G (7), ?/l(7)]} 
7 e [1,00), <;5* (7) < 771,(7) <7;|j 



^1(7) = ^ 1/2 ^ I' (4.3.4) 

where C,, (to be determined in Proposition 14 . 1 may depend on 7 but not on 



Moreover we restrict 77* (7) so that (|1.3.ip is satisfied. 

Lemma 4.3.2. For sufficiently small e, with bound depending on Cn'- 

There is an 7?£(7) G [0:'7l(7))j continuous in 7, so that hi .3.1]) is satisfied for 
all Vl e ivlil),V*n}' l>'^i never for tjI G (0,?7£(7)]. 
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Figure 14: yy^ below 77^ violate (|1.3.ip . For 7 > 1 some distance from rj'^ is 
needed for technical reasons. The Leray-Schauder degree argument uses param- 
eters in the dark shaded set A which is path-connected for small e > 0. 

Proof. (|1.3.ip is satisfied for ij^ ~ rj*^ (see Figure IH]), by Proposition 12.6.1] in 
this case L shock and R shock coincide, and the R shock never intersects the 
circle with center = vi and radius c„ = c/. 

Clearly the distance to that circle is strictly decreasing as ry£ decreases. By 
continuity there must be an rjf^, depending continuously on 7, so that for 77^ = ry^ 
the L shock touches the circle. For all smaller 77^ > the circle is intersected, 
because JVij^ylj > does not allow the L shock to pass below the circle. 

Since the L shock tangent and location depends continuously on 7, ry£ must also 
be continuous in it. □ 

Lemma 4.3.3. For sufficiently small e, with hound depending on Cn'- 
For all rj*^ (7) ( continuous in ^) that satisfy 

vUj)<vlil)<vlh), (4.3.5) 
A is path- connected and contains the point (7,772) = (1,77^). 

Proof. The curve 7 1— *■ ■^ti*j^ + tj*l) is continuous and contained in A. It intersects 
every line {7} x [jy^i^Lli including the one for 7 = 1 which is (77* , 77]^]. A is the 
union of all these lines, so the proof is complete. □ 

4.4 Function set and iteration 

Definition 4.4.1. Let U C IR" open nonempty bounded with dU uniformly 
Lipschitz. Let F C dU. For k € INq, a G [0, 1] and /3 S (—00, k + a] we define 
the weighted Holder space Cp'"{U, F) as the set of m e C'='"(J7 - F) so that 
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is finite. 

Definition 4.4.2. For sufRcieiitly small 5p > 0, there is a function b e C^(IR^) 

with 6, |V&| < 1 so that 6 = on and & > elsewhere, b even m 77. 

From now on we fix a particular b. 

Proof. The construction is straightforward, dp is taken so small that Pj^^ and 
P^^ are well-separated. That is possible because ^2^°^ < ^^^^ . □ 

Definition 4.4.3. 
Onion coordinates 

To define a function subset in a fixed Banach space, we need to map the do- 
main with its free shock boundary to a fixed square. Define a C°° change to 
coordinates (cr, r;) S IR^ (see Figure [T5|) so that 

1. 77 is preserved and (T = (T{£,,ri) aC°° function with a function^ = £,((J,ri) 
inverting it, 

2. P^'' maps to a subset of {a = 0}, 

3. Pjj"^'' maps to a subset of {a ~ 1}, 

4. 41 maps to (0,0), 

5. maps to (1,0), 
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6. and A^'^^ maps to {a G (0, 1), 77 = 0} precisely. 

7. Let T C 5^ be the set of unit tangents of P^^ and P^^'' ■ For some constant 
dp„ > depending only and continuously on A G A, require that for every 
unit tangent t of a cr level set, 

d{t,T)<Sp,. (4.4.1) 



We require that the change of coordinates depends continuously (in C°°) on 
A S A. The construction is straightforward. 

Here and in what follows, we will use the weighted Holder spaces C|'"(C/), as 
in Definition HXlJ The domain U is either [0, 1]^ with F = {(0, 1), (1, 1)}, or 
n with F = {£,L,^B.} (to be defined). For the shock we use U — [0, f ] with 
F = {0,1}, or U = [^L,Cfl] with F = {^l,£,r}', for the arcs only the upper 
endpoints are in F and for the wall A we have F = %. We omit F as it will 
be clear from the context. C^'" are Banach spaces so that standard functional 

analysis applies. Moreover, C|'"(rj) is continuously embedded in C^(0), so 
we have regularity in the corners as well, which is crucial. (3 £ (1,2) and 
a € (0, — 1] will be determined later. 

Free boundary fit 

Let T be the set of functions -0 € C|'"([0, 1]^) that satisfy all of the many 
conditions explained below. Require 

U\\cr^[^,i]^)<Cc{e). (4.4.2) 



For all (T e [0, 1] define 

V>(a,l)-V>^(0,0) 
s cr) := y ; 4.4.3 

f J 

it satisfies ^{(t,1) = ?/;^(^(cr, 77), 77) with 77 = s(cr). We define another coordinate 
transform by first mapping (cr, C) G [0,1] to ((7,77) with 77 = C*('''): ^i^d then 
mapping to ^ with the previous coordinate transform. 

Let resp. be the ^ coordinates for the (cr, C,) plane points (0, 1) resp. (1, 1). 
Let S be the ^ plane curve for (0, 1) x {1}. Define Pl, to be the images of 
{0} X (0, 1) resp. {1} x (0, 1). Finally, let VL be the image of (0, 1) x (0, 1). 

Require shock-wall separation: 

max s(cr) > 8sA > 0. (4.4.4) 

<t6[0,1] 

Require 

^ — ^(cr, C) is strictly increasing in cr for C = 1. (4.4.5) 
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Require: corners close to target: 

la-ClU6^-G,l<el/^ (4.4.6) 

We require e to be so small that G Pc {C — L, R), i.e. may not be higher 
than the upper endpoint of Pc- 

For later use we define := r ^ ^^^^ ''^^'^ ^"^^ ?l b. ^^^^ ^lr ^ 

Corner cones: require 

sup_Z(e-Cc,^'-Sc)<7r-fcc (4.4.7) 

(As discussed in the introduction, Z{x,y) is the counterclockwise angle from x 
to y.) 

(|4.4.4p ensures that the map from (cr, C) to ^ is a well-defined change of coordi- 
nates, uniformly nondegenerate (depending on 6sa and Cc), with C^'"([0, 1]^) 
resp. C^'"(ri) regularity. Since the step from {a,() to (cr, ?y) uses ip, the entire 
coordinate change is as smooth as V'- If we prove higher regularity for ip either 
in (cr, C) or (^, ry) coordinates, we immediately obtain the same higher regularity 
for the coordinate transform and in the respective other coordinates. 

It is clear now that dil is the union of the disjoint sets S, Pl, Pr, A, {£,r}, 
{^l}, {^bl} and {^b_r}. By (|4.4.4p . is a simply connected set. 

(|4.4.5|) ensures that s can be defined as a function of 5, which is the way we use 
it from now on. 

Iteration 

Shock strength/density: require that 

-X-^|Vx|'>0, (4.4.8) 

so that p is well-defined, and require 

min p > pi + Sp. (4.4.9) 
n 

Pseudo-Mach number bound: require 

<l-SLb-b in n, (4.4.10) 

(Note that L is well-defined because by (|4.4.9p p > 0, so c> 0.) 5 = on P|°)j 
which have distance > | (for sufficiently small e) from 51, so (|4.4.10p implies 

<l-l\Vb\L^SLb-e<l-l6Lb€ inH, (4.4.11) 
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Require: there a function -0 € C^'"(f7) with the following properties: 

1. ip close to V': 

where r/ G C(JF; (0,oo)) is a continuous function to be determined later. 
Here and later we regard ip as defined on [0,1]^ instead of fl, via the 
coordinate transform from ^ to (cr, C) defined by "0 (see above) . 

2. We require r/ to be so small that 

-X-^|Vxl'>0, (4.4.13) 

V^Jj^vi, (4.4.14) 
so that in particular p is well-defined and positive. 

3. Moreover we require 77 to be so small that 

{4 + (1 - 7)(X + ^iVxH)/ - Vx' > 0, (4.4.15) 
i.e. is a (symmetric) positive definite matrix. 

4. Let C = £{ip,ip) be defined as 

((c2 + (1 - 7)(x + ilVxH)/ - Vx") : V^V', (4.4.16) 
iVxP , (l-e)((7-l)x + c^) 



2 2+(l-e)(7-l) ' 



(4.4.17) 



(pVx - P/VxO • ^ti ' (4.4.18) 

0,). (4.4.19) 

(|4.4.18p is well-defined by (|4.4.13p and (|4.4.14p . The other components 
have no singularities. 

e so |VxP e C^L"i, so ((c2 + (1 - 7)(x + i| VxP))/ - Vx^) S 

C^jfi CO'/'-i C"-" (note a < - 1 as required above), and V'^ip S 
C;^-"2. so gXTO is e C;^-:^2. In the same way we check that (j4.4.17p . 
(|4.4.1g|) and ()4.4.19P are C^'"^. Hence we may take the codomain of C to 
be the Banach space 



^^ip is the product of an iteration step with input ^. We will show/ensure in Proposition 
I4.4.7l that ip is unique and continuously dependent on ip. 
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Alternatively, if we consider the puUback to (cr, C) coordinates, as defined 
by ip above, we may consider 

C°L"2([0, 1]') X C^L"i [0, 1] X 4_"i[0, 1] X Cl-^^,[0, 1] x 4_"J0, 1]. 

In the same way wc can discuss ip cither in C|'"(ri) or in C^'"([0, 1]^). 

With these topologies, clearly £ is a smooth function of ip and i/j. 
Most importantly: require 

£(V',i/') = 0. (4.4.20) 

Other bounds 

Require 

||V^|lco,i(n) < Cl (4.4.21) 
where may not depend on e. 
Xt and Xn on parabolic arc: 

JB^2^c''\% < Cpte'/^ (4.4.22) 

max c-i^ < -Spn. (4.4.23) 
PlUPr on 

Here, 6pt, Spn may depend only on A, but not on e (or t/j). 
Horizontal velocity: 

maxi/'c < (^1,^61/2. (4.4.24) 
fi 

Vertical velocity: 

V'r, > + Syy in n. (4.4.25) 
Left corner shock tangential velocity: 

rnax Vt/i x Hl > vl >^ nL + C^Lf^l"^ (4.4.26) 

Shock normal: Let N C (unit circle) be the set from Tin to ul (counterclock- 
wise) . 

sup d(n, TV) < Cs„e^/^ (4.4.27) 
s 
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Set I]i := Pl, ^2 := S, := Pr, and E4 := A. Write the components (|4.4.17p . 
(14.4.191) . (14.4.181) of/: as 

3XC,x(0,Vx(0) (* = 1,...,4), 



where the ^ dependence includes the dependence on x(^) and Vx(C)- 

has some singularities, but not on the set of ^, x, Vx so that ()4.4.25p and 
(|4X9l) (resp. ()4.4.13p and (j4.4.14p ) are satisfied. That set is simply connected, 
so we can modify on its complement and extend it smoothly to O x IR x IR^. 
The modification is chosen to depend smoothly on A. 

Require uniform obliquenesJ^: 

\9yn\>So\g},\ VeeS\ (4.4.28) 



Functional independence in upper corners: for i,j = 1, 2 and for i,j = 2, 3 set 

J 



G 



regard it as a function of ^ (including the dependence on Vx(0) ^^'^ require 
\\Gl\\G-^\\<Cd in BsA^c) nn, G ^ L, R. (4.4.29) 

Let T be the selfH of admissible functions so that all of these conditions are 
satisfied. Define J- to be the set of admissible functions such that all of these 
conditions are satisfied with strict inequalities, i.e. replace <,> by <,>, "in- 
creasing" by "strictly increasing" etc. 

[This is the end of Definition gXl] 

Remark AAA. If -0 = ■(/;, then (|4.4.16p . (|4.4.18p . (|4.4.17p . (|4.4.19p and the defi- 
nition of S yield 

{c^I - Vx^) : VV = inU, 

X,, = on A, 

X^ = X and 

(pVx - P/Vx^) • n = on ^, 



L = ^fY~e on Pi U Pr 
(we may take closures by regularity (|4.4.2p '). 



^^In various articles Lieberman uses the term "uniformly oblique", probably with "oblique" 
in the sense of "not tangential" . We adopt this terminology instead of the more common but 
less useful sense "not normal" (see e.g. IGT83I or [PP97])- 

^^The notation T docs not necessarily imply that T is the closure of T . 
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Remark 4.4.5. In any point on A, we can use even reflection (see Remark l2.2.ip 
of "0 across A to obtain a point in the interior, or (in the bottom corners) a point 
at a quasi-parabohc arc with the interior equation applying inside. on A, 

for even reflection of ip, implies that the solution is across A; then necessarily 
it is also C^'" (away from the shock-arc corners). 

For ip standard regularity theory immediately yields that the solution is analytic 
in the bigger domain near A. The same technique applied to ■0 and to solutions 
ip of linearized equations (here ip, ip and tp are reflected) yields C^'" regularity 
(away from the shock-arc corners). 

Proposition 4.4.6. For sufficiently small e (with bound depending only on Cpt ) 
and rj (depending continuously and only on ip^S^y): 

for all ip & T, £(V',V'') well-defined for ij/ nearip, and the Frechet derivative 
dC/ dip)' {ip,ip) (of C with respect to its second argument ip' , evaluated atip' = ip) 
is a linear isomorphism of C^'"' onto Y. 

Proof (|4.4.18|) is the only part of C with a singularity. (|4.4.9p and (|4.4.14p 
guarantee that it is well-defined and smooth for ip = ip. For -0 in a sufficiently 
small neighbourhood of ip (i.e. take r/ small), it stays well-defined. 

Let tp € C^'"(r2); it is meant to be the first variation of 'ip (and, at the same 

time, x). dC/dip{ip,'4>)'ip is a tuple of fmictions in Y. Its first component (see 
(|4.4.16p ) is of type 

(c^/- Vx^) : V^V^-f-fo- VV^ (4.4.30) 
where b is some vector field. 

We linearize (|4.4.18p : here we can use that we linearize a.t ip = ip, so 

\vi - Vip\~'^{vi - Vip) = n 
because ip = ip^ defines the shock (see (|4.4.3p in Definition 14. 4. 3p . Result: 
p{l - c~^xl)-pn - pc'^XnXfipt - pc'^ip 

- [pVx - Pi^x'Wi - VVr'(l - (^|— ^)')VV^. (4.4.31) 

—n 

The coefficient p{l - c~^Xn) of V'n: positive by (|4.4.9p and ()4.4.1ip . has the 
opposite sign to the nonzero coefficient — pc^^ of ip, negative by (j4.4.9p . 

The other components linearize to 

Vx • VVa (parabolic) (4.4.32) 
X„. (wall) (4.4.33) 
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Note that the coefficient vectors of Vip are the same as in (|4.4.28p and 
(|4X29ll . 

To prove the Proposition, we apply [Lie88b[ Theorem 1.4]. Our C^'"(r2;i^) 

spaces correspond to his spaces h'^J^ if a > and F C dO, is finite and 
contains all points where dO. is not — as is the case here. We check the 
preconditions: 

1. Yii are curves (except perhaps for the endpoints), meeting in single 
points (corners) . In each corner the two curves meet at an angle < 9ij < 
Of), with 6q < tt. (0 < % is trivial from (|4.4.27p ). The lower bound is 
obvious; the upper bound is obvious for corners with A and supplied by 
(|4.4.7p for corners with S. 

2. The equation is uniformly elliptic, by (|4.4.1ip . 

3. All boundary operators arc uniformly oblique, by (|4.4.28[) . 

4. Condition (1.17) in loc.cit. with /3i ^ (3j is equivalent to (|4.4.29p . 

5. Finally, the only appearance of tp is in (|4.4.3ip . where its coefficient his 
nonzero and has the opposite sign as the coefficient of tpn (note that n 
is the inward normal here). (This observation allows to apply maximum 
principles to obtain uniqueness.) This means (1.19) in loc.cit. is satisfied. 

6. The bottom corners arc not a concern: by Remark 1 4 . 4 . 5 1 we can use even 
reflection across A. 

7. The other preconditions are technical and easy to verify. 

Theorem 1.4 in loc.cit. yields that dC/dx is an isomorphism on C^'" onto Y, if 
we choose a £ (0, 1) and /? e (1, 2) sufficiently small, depending on the constants 
So, Cd, 5Lb: 5d and Scc- □ 

Proposition 4.4.7. r/ can he chosen so that ip is unique and depends contin- 
uously on X and on ip G (both ip ad ip in the C^'"([0, 1]^) topology). 

Proof. We use the subscript A for C, T here to indicate their dependence on it. 

Take r/ := 1 first. £a(V'iV'') is well-defined for all '\\} E T\ and t/"' = ^A, as 
well as sufficiently small perturbations of , by (|4.4.8p and (|4.4.1ip . It is easy 
to check that there is an r2 > 0, depending continuously on -0 and A, so that 
ii' £ Br^ii^) ^ C\{^,ijj') is well-defined and . Take r/ ^ min{r/,r2}. This 
may shrink !F\, but the properties of C are not changed. 

Consider a particular A G A and a corresponding ip G T\. By Proposition 
14.4.61 and the inverse/implicit function theorem for Banach spaces, there is an 
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r e iO,ri] so that if,'' e -Br(V') ^vii^'ji^') is a difFeomorphism for every 
tp' e 5^(7/-) and A' e A n Br{X). 

Let ra be the supremum of all r with this property. (-0, A) 1— > r^^ip^ A) is contin- 
uous: for any V'" G S^IV-) and A" G Br(A), set r" := r-max{|V'-i/'"|, |A- A"|}. 
Then Br"{^") C and B,-"(A") C Br(A), so V^' € Br-{ip") ^ Cx'{ij',4'') 

is a diffeomorphism for aU e Br"{ip") and A' e A n Br"{X"). Therefore 

r3(V'", A") > r3(^, A) - max{|^ - |A - A"|}. 

On the other hand, if max{\4) - iP"\, |A - A"|} < ir3(i/', A), then r3(V'", A") > 
r3(?/', A)/2, so V' G Br-3W,",A")(V'") and A G Br^(^^„^x"){>^"), so we may apply the 
same argument with roles reversed to obtain 

rsi^, A) > r^ii;", A") - max{|^ - ^"|, |A - A"|}. 

Clearly is continuous. 

For r = r3, the property need not hold, but we take rj -f— min{r/, ^r^}. This 
may shrink ^ more, but again the properties of £ and r2 above are not changed. 
With this choice, ip from Definition l4.4.3l must be unique (determined by (|4.4.12p 
and (|4.4.20p ). It is also clear from the properties above that (V', A) 1-^ i/' is a 
continuous (in fact C^) map. □ 

Proposition 4.4.8. For Sti, e rj sufficiently small: For all continuous 
paths t G [0,1] i-^ X{t) in A, Utg(o,i) ^ ^Ht)) °P^^ Ute[o,i] ({0 x 
Tx{t)) is closed in [0, 1] x C^-"([0, l]^). 

Proof. All conditions on ■0 in Definition 14.4.31 are inequalities which can be 
made scalar by taking a suitable supremum or infimum. Then their sides are 
continuous under C^'"([0, 1]^) changes to -0, and therefore -0. (Most inequalities 
need only Ci([0, 1]^).) 

1. Closedncss: consider sequences (tmipn) in Ute[o 1] (i^i ^ -^K*-)) ^^^^ con- 
verge to a limit (t, tp). 

Let tpn be associated to ■0„ as in Definition 14.4.31 By continuity (Propo- 
sition [523), ii^n) converges to a limit tp as well. By continuity oi C in ip, 
ip and A, we have C\(t){ip, = as well. 

By (|4.4.3|) for s„, i/)„ instead of s, ip, (s„) converges in C^'"[0, 1] as well, to 
a limit ,s which satisfies (|4.4.3p itself. 

^^We make no statement about T being the closure of T. It certainly contains the closure, 
but it could be bigger, for example if one of the inequalities in Definition 14.4.31 becomes 
nonstrict in the interior without being violated. 
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Most conditions on 'ip are nonstrict and continuous inequalities, so they 
are still satisfied by ip. We check the strict inequalities explicitly and in 
order: 

(|4.4.5p : this is implied by (|4.4.27p and the nondegeneracy of the {a, () 
77) coordinate change. 

(|4X8ll is imphed by ([iXg]) . 

(|4.4.13p resp. (|4.4.14p resp. (|4.4.15p are implied by (|4.4.12p resp. (|4.4.25p 
resp. (|4.4.1ip . by choosing rj sufficiently small. 

All inequalities arc satisfied, so G !F. 

2. Openness: same proof, using that all inequalities are strict now, by defi- 
nition of JF, hence preserved by sufficiently small perturbations. 



□ 

Definition 4.4. 9. Define K. : T ^ C^'"([0,l]2) to map ^ into V' as given in 
Definition l4.4.3i but pulled back to {a, Q coordinates and the [0, 1]^ domain (see 
Definition 14. 4. 3p with the coordinate transform defined by ip- 



4.5 Regularity and compactness 

To obtain regularity at the free boundary, wc need a kludge: a transformation 
to a fixed boundary problem. With some further advances in elliptic theory this 
step should become obsolete. 

Remark 4.5.1. Using the ((7,77) coordinates from Definition 14.4.31 define the 
coordinate transformation from ^ to x = (cr, 77) and then to y ~ (a, fi) where 

fj, := i^{cr, rj) — tj/ {a, rf) = -0(0-, rf) — 4'^ {Q, 0) — v^j-q. 

ifj = ^p^ on the shock, so = there; fi > in ft (by (I4.4.25|l ). 

We check that this transformation satisfies all conditions of Proposition 15.2.11 
and Remark 15. 2. 2| where we take u = ip. The ^ to a; part is trivial, since it is 
independent of 4. We study the second part: 



Vx 



Vu = 



9-0 



= [V'ct tp,,] 



Obviously yx is uniformly invertible, and 



VuU^ 



1 
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is uniformly invertible by (|4.4.25p . The norms are bounded by constants de- 
pending only on 5r,y (to bound i/'?) ~ below away from zero) and Cl (to 
bound ij}, V-0). 

All conditions of Proposition 15 . 2 . Il and Remark 15.2.21 are satisfied, so we obtain 
the analogues of (|4.4.1ip . (|4.4.28p and (|4.4.29p : the coordinate transformation 
yields a uniformly elliptic equation with uniformly oblique boundary condi- 
tions that arc uniformly functionally independent (as functions of Vip) near the 
corner. The constant for each property grows at most by a factor depending 
continuously and only on S^jy ,Cl- 

Proposition 4.5.2. For sufficiently small [3 G (1,2) and a € {0,(3— 1), de- 
pending only on Cd^Shb ■ e, <5o, Cl, 5vy: 

1. 

\\4co.^[U,U] < CsL (4.5.1) 



\\4cria,U\^C. (4.5.2) 

for CsL = CsLiCL,Svy) and Cs = Cs{Cc,Svy). 

2. For a fixed point of K,: 

(a) \4.4.21^ is strict for sufficiently large Cl- 

(b) l4-4.2^ is strict for sufficiently large Cc = Cc{Gd,5Lb'^,GL,5o,5vy,^d)- 

(c) For K ^ f7 - {11,6?} and all k > 0, a' E (0, 1), 

Mc^..'^K)<CcK (4.5.3) 

where Cqk = CcKid,CL,So,Svy) is decreasing in d := d{K,P^^ U 
P^^) and not dependent on e. 

(d) ijj is analytic in fl — {S.l,^r}- 

3. For sufficiently small rj > 0, depending continuously and only on ip, there 
are (5a, (5/3 > so that for all i/j E 

U^..s.^^^<C^ (4.5.4) 
Here, C]c,Sa,5p depend only on Cd, Shb ■ e, (5o, CL-,5^y, 
Proof. 1. The shock is implicitly defined by 

V'(c, 5(c)) -v/(e, 5(0) -0 (ee KL,ei?]) 
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(compare Definition 14.4.31 (|4.4.3p '). The derivative with respect to s of 
the left-hand side is V'r; — Wj. (|4.4.25p bounds this away from zero. The 
implicit function theorem yields the derivative part of (|4.5.ip : bounds on s 
itself are supplied by integration and the part of (|4.4.6|) (for example) . 
The implicit function theorem also yields (|4.5.2p from (|4.4.2p . Henceforth 
we use without further mention that estimates on x automatically yield 
corresponding estimates on s. 

2. Now consider a fixed point i/' = /C («/>). 

(a) (|4.4.1ip implies L < 1 which can be rewritten 

IVxP < M±(i_2k). ,4.5.5) 

7+1 

xi^B.) — X^{(,r) which is bounded by (|4.4.6p . Using the differential 
inequality (|4.5.5p vertically downwards from to A, then along A, 
and finally along any vertical line upwards to 5, wc achieve uniform 
bounds on x S'lid |Vx| in all of fi. This is (|4.4.2ip which is strict if 
we take Cl large enough. 

(b) By (14. 4. lip the equation 

(c^/- Vx^) : = 

is uniformly elliptic. Standard dc Giorgi-Nash-Moscr and then Schauder 
theory converts (|4.4.2ip into interior C'^'" estimates for any k > 2 and 
a e (0,1). 

At the shock and near the shock-parabolic corners we first apply the 
coordinate transform from Remark 14.5.11 This is necessary because 
we have only a C°'^ bound of S at this point; using C^'^ regularity 
would cause circularity. The transform yields a new problem with 
fixed boundary (both arc and shock are mapped into straight line 
segments). The uniform cllipticity (|4.4.1ip . the uniform obliqueness 
of the boundary conditions (|4.4.28p . and their uniform functional 
independence in the corner (|4.4.29p are preserved by the transform. 
It is also crucial that the boundary conditions themselves are smootli. 
For this purpose we have combinetQ the two shock conditions (|2.3.ip 
and (|2.3.2p so that the boundary operator in (|4.4.18p is smooth in ^, 
x(^) and Vx(^). This also requires a bound of p below away from 
zero, by (|4.4.9p . as well as a bound of away from vj, by (|4.4.25p . 
At the boundaries away from the corners, we use [Lie84| . which yields 
(j2,a regularity for any a £ (0, 1). Near the corners we apply Propo- 
sition [5TTTT] It yields C|'" regularity for some range of a S (0, 1) and 
(3 e (1,2); we fix a,/3 below. 

'^■'If we had used 12.3.21 1 by itself, it would contain n for which we have only an L°° bound 
at this point. 
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We obtain (|4.4.2p . with Cc = CcjCd, SLb-e, Cl,So, 5.vy,5d)- 5yy is due 
to the transform from Rcmark l4.5.1l Moreover Cs = Cs{Cc,5vy), as 
discussed above. 

(c) In every K (s^ Vl with positive distance to Pl and Pr, we use that 
(|4.4.10p provides 

sup LF' < 1 — Slij sup b 

K K 

where the right-hand side is < 1 and independent of e (by Definition 
I4.4.2p . Thus the equation is uniformly eUiptic in K, with eUipticity 

constant depending on d{K^ pf^ U P^ii), but independent of e. Since 
the obliqueness constant So in (|4.4.28p was already independent of e, 
the same arguments as before yield uniform in e regularity in each 
K, for any Holder norm. This is (|4.5.3p . 

(d) For elliptic problems with analytic cofhcients, analyticity of the so- 
lution is classical (e.g. [Mor66[ Theorem 6.7.6']) in the case of fixed 
boundaries. To deal with a free boundary, we use the transform from 
Remark 14.5.11 again . Loc.cit. yields analyticity for the new problem. 
The inverse coordinate transformation is defined in terms of the new 
coordinates and the solution of the new problem, so it is analytic as 
well. Then tp itself is analytic in il — {£,l,^r}- 

3. For a general ip ^ !F, not necessarily a fixed point of /C, and ip h^ ^ are 
different. The boundary conditions arc oblique derivative^ conditions: 

Each g'' is C^'" in f (fron{lf| x in (|4.4.17|) ). in particular C^^'^'^ 

Moreover each g'^ is C°° m^{£^) and V'0(^). For (|4.4.18p this requires that 
p is uniformly bounded below away from 0, and Vijj uniformly bounded 
away from vj. This is guaranteed by (|4.4.9p and (|4.4.25p (bounding p rcsp. 
VV') combined with (|4.4.12p (bounding \p — p\ resp. |VV' — VV"!), if rj is 
chosen small enough (depending only and continuously on 'ijj and 6yy). 

The boundarief0 P'^ are A,Pl,Pr^ which are perfectly smooth, and S, 
has a C^'" bound (l4X2l) . implying C^^P . 

Uniform eUipticity, obliqueness and functional independence still hold by 
(|4.4.1ip . (|4.4.28p and (|4.4.29p . combined with (|4.4.12p . for r/ sufficiently 
small(er). 

All combined, Proposition 15.1.1) and [Lie88a| Corollary 1.4] (away from 
the corners) yield C^'^ regularity. Here, k depends on /3 only whereas A 

For a Dirichlet boundary condition, C^'^ data would yield C^'^ , but not C^-l^+^fi regu- 
larity. But for a derivative condition we can gain one order. 

^^Here it is crucial that V'i/' does not appear anywhere in 114.4. 171 1. I l4.4.18t . 14.4.191 1: other- 
wise 3* would only be C^'" in ^ which is not enough to gain regularity. 

^^Here we can apply Proposition lS. 1 .Tl without the coordinate transform from Remark l4.5.1l 
because the shock is already known to be C^'''. 
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does not depend on a, j3 at all. Therefore we may pick /? = 1 + A/2 and 
then a = k/2 for the resulting k. With Sp = A/2 and 5^ ~ n/2, (|4.5.4p is 
satisfied for suitable C]c. 

□ 

Proposition 4.5.3. For sufficiently small r/ > (depending continuously and 
only on tp, S^y), K, is a continuous and compact function of ip and A on IJ;^^^ 

Proof. Continuity: we have already shown in Proposition 14.4.61 that -0 is a 
continuous function of ip. Therefore /C depends continuously on A and ?/;. 

Moreover JC is compact: by (|4.5.4p in Proposition 14.5.21 the range of ^/i is a 
bounded subset of (fi) which is pre-compact in C|'"(f7). Fullback to a, C, 

coordinates is continuous in the latter topology, so the image under it is still 
pre-compact. Altogether /C is a compact map. □ 



4.6 Pseudo-Mach number control 

Proposition 4.6.1. For e and S^h sufficiently small, with bounds depending 
only on 6p: if tp Cz J- is a fixed point of K,, then \^.4-.10^ is strict and 

L^<l-e inTi-PL-PR. (4.6.1) 

Proof. 

d(n,pf upr)>l.e, 

for e small enough. Remember from Definition gXH that & = on pf U P^^ . 
Therefore, on Pl U Pr: 

L^ = l-e<l- \\b\\co.i ■ d{PL U Pr, PP UP^°'^)<1- 5Lb ■ b, 
<i 

e.g. for 6Lb <1. 

On the shock, we may use (|4.4.9p combined with Proposition 13.6.11 to rule out 
that + (Sif, • 6 has a maximum in a point where L < 1 and L > 1 — Sls- 
Here SLb has to be chosen so that \SLb^b\ < 6ls is satisfied. (Now S^b depends 
continuously on Sp as well.) Then + SLb • 6 < 1 as well if, again, S^b is small 
enough. 

In addition we can choose dLb so small that Shb ■ b satisfies the preconditions 
of Theorem 1 and Theorem 2 in [EL05aj (where it is called 6). Note that b 

^*We may actually lose regularity because the puUback is defined by tp, hence only C^'", 
but we do retain compactness. 
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Possible values 
in 6? 



Figure 16: x, Vx on Pr define a patli (p {p{4')j ^(0)) (bold curve). The path 
is constrained by (|4.7.1ip and (|4.7.10p . The possible k,p in are on the thin 
solid curve which is parametrized by r/j^; one endpoint of the path is on this 
curve. The other endpoint is on the k axis. 

is even in rj by Definition I4.4.2[ so db/dn = on the wall. Let Slq be the d 
from those theorems (it depends only and continuously on A). Then +SLb • b 
cannot have a maximum in a point where > 1 — S^u- If, again, Sib is chosen 
sufficiently small (no new dependencies), then + 6^, • 6 < 1 in U A, hence 
in n. Therefore (|4.4.10p is strict. 

(|4.6.1[) can be shown in the same manner, by taking 6 = instead, using the 
actual boundary condition L — VT — e on Pl , Pr and and considering e < 
SlStSlu- □ 

4.7 Second derivatives on arcs 

Consider a fixed point x G of /C, so that = 1 — e is satisfied on Pl,Pr. 
This can be restated as 



In what follows, ji will be an outer normal to P^, or Pr and t — n^; i9„ and dt 
are defined accordingly. When considering a particular point on an arc we also 
use r and v, which are Cartesian coordinates such that dr = dt and di, = dn in 
that point (only). 




(4.7.1) 
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Take dt of (HTID (using dtdr = and dtdy = drdy): 

(7-l)(l-e) 
= XtXtt + Xi^Xiyr H — -, -tXt 

7 + 1 -£(7-1) 

2 

= XrV'rr + Xi/V'i^r — -, "tXt- 

7 + 1-6(7-1) 



(4.7.2) 



We have two equations ( (14.7.2^ and the interior equation (|2.2.4p ) for the three 
components of -D^X- That is the case for other equation and boundary con- 
ditions as well; it is not really sufficient for any control. But here there is an 
additional tool: (|4.6.ip implies that L attains its maximum in f2 in every point 
on PlU Pr. Therefore {L'^)u > on Pj, U Pr, or equivalcntly 



7 + l-e(7-l) 



Xu>0. 



(4.7.3) 



From here on — and in the next few sections — we focus on Pr; all results have 
an analogous extension to P^. 

Combine gZ31), (SZll), (PXSD : 



Xr 





Xi' 

Xr 





Xi^ 



Xr -2XtXu - xi. 



7 + 1-6(7-1) 



Xr 




"0' 


Xu 


+ 


a 











(4.7.4) 

where we use a "slack variable" a > 0. Solutions (we use x^ = (1 ^ e)c^ ~ X?)- 

(4.7.5) 



tprr = 



2 n + e xl \ (c'-x')x. 



7 + 1 — e(7 — 1) V 1 — e c" 



(l-6)c4 



The a term is positive due to Xn < (by (|4.4.23p ). 

Now consider radial coordinates, centered in the origin, with = correspond- 
ing to the positive ^ axis. Pr is at a fixed radius r = (1 — e)^/^Cfl and covers 
4) G [0, (/>]. In these coordinates: 

X00 = r'^iXrr - r~^Xy) = r^ii^rr - 1 - r^^Xiy) > f{c^,X^), 

2 .1 + e p2 



7+l-e(7-l)U-e h 



- er^) - + ^Jr'^h{\ - e) - 



On the other hand, = + (1 - 7)(x + |VxP/2) and (|4.7.ip yield 



(c^)./- = 5(c^,X0), g{h,p) 



-2(7-1) 
7 + 1-6(7-1)^ 



(4.7.6) 



(4.7.7) 
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We seek stationary points of the ODE system {p^,h^) = {f,g)- g{ho,po) = 
obviously requires po = 0. 

2e 



O^f{ho,0) = -r' 1 + ^ 7 TT +rVhoVl 



7 + l-e(7-l) 

n \ 2 2 



e 



^ ho= [1 + -, -T (4.7.8) 

g is already linear; we linearize / around the stationary point: 



dhfiho, 0) = r J'^ = , — Te T ^- 



2y/ho 2(1 



7 



2e 

+ l-c(7-l)' 



For easier treatment of the isentropic case, we change coordinates again (see 
Figure [TB|) : take 

k:^{-aj/agy^\h-ho). (4.7.9) 
From linearization we have 
P0 > <Jfih - ho) + 0{{h - hof+p^) = {-(jfagY'^^ -k + 0(fc2 +p^)^ (4.7.10) 

= :ae 

h = {-(^f/^gY'^^aP = -~{'^f(^9f'^P = -(^ep. (4.7.11) 

Here and later, 0(a;") is a term with absolute value < Cx", for \x\ < R, where 
C and R may depend only on Cl, but none of the other constants. While the 
elliptic equation degenerates as e J, 0, all ODE we discuss here are well-behaved 
for e I 0. However, some 7 i 1 are dehcate, because ag = 7 — 1 + 0(e), so (|4.7.9p 
has a singularity; we make detailed comments in each case. 

Represent {p, k) by radial coordinates (g, 6) , with 9 = corresponding to the 
positive p axis and 9 = ^ to the positive k axis. Then 



pk^ - kp^ _ j < 0{q) - gfl, fc > 0, 
g2 ^\>0{q)-ae, k < 0. 



(4.7.13) 



4.8 Arc control 



Proposition 4.8.1. If Cpt < 00 is sufficiently large, if Spn > is sufficiently 
small, if e is sufficiently small and Cpy , Cpp sufficiently large, with bounds de- 
pending on Cpt, then for any fixed point x ofIC, ^J^JTS^ and ^J^J72^ are strict. 
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and 

\p- Pc\<Cppe^''^ and (4.8.1) 
\v-vc\<Cp,^'^ onPc (C = L,R), (4.8.2) 
and in the case 7 > 1; 

?(f,Y-^9^) + 2^2 onPR, (4.8.3) 



(-| + Te0, f) + 27rZ onPL 



Proof. We focus on Pp first. 

For ~ we have p{(j)) — and c{(f)) ~ cp, hy construction of the R shock 
in Section p{(j>),c{(p) depend smoothly on ^p, so (|4.4.6p yields 

\pm < (4.8.4) 
|c(^) - cr\ < (4.8.5) 

in ^/f, for some constant Cpk- = 1 — e, combined with (|4.4.22p . as well as 
(|4.4.22p integrated in tangential direction, then implies 

?((/))< C,ei/2 V0e[O>]. (4.8.6) 

for Cq ~ Cq{Cpt). (C"^ may he > Cpt, so this does not imply the sharp form of 
(liX^ yet.) 



Isothermal case For 7 = 1 we need to consider only a single differential 
inequality, since c is constant. (|4.7.6p takes the form 

P0 > -(l-e)c' -6 + 0(^2). (4.8.7) 

^ = is the corner between A, Pp where x<t> = fXri = by the boundary 
condition on A, so p{4> = 0) = 0. Integrating (|4.8.7p from to yields 

p(0) > 0{e). (4.8.8) 

On the other hand p(</)) is controlled by (|4.8.4p . Integrating (|4.8.7p from (j) to (p 
yields 

p{(t>) < 0(e^/2). (4.8.9) 

(|4.8.8p and (|4.8.9p combine to 

c"V"^ max p{(j)) < Cpte^/'^ 

0G[O,c#.] 

for e sufficiently small and Cpt sufficiently large. This is the strict form of 
(|4.4.22|1 . 
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Isentropic case For 7 > 1: first we show (|4.8.3p . We fix 6* £ [0, 27r) here. 

Assume that q{<j)') > and 9{<p') £ (tt, ^ — agcj)) for some 0' G [0, 0]. Necessarily 
> because p(0 = 0) = 0, so either 0) = or = 0) = tt ± f . 

Let 00 S [0, 0') be maximal so that 



^'0 



)^(^,^) or g(<^o) = 



Such a 00 must exist because p(0 = 0) = 0. For e (0o,0'), fc(0) < and 
pW < 0. 

By the p < part of (|4.7.12p (in reverse direction), q{(p) > for e [0o,0']- 
On the other hand (|4.8.6p applies. Therefore q{4>o) = is not possible, so either 

eiM^TTOV 0(0o) = f . 

(|4.7.1ip with p < (in reverse direction) shows that k{(j)o) < k{<j>') < 0, so 
9{(j>o) = TT is not possible. 

For e sufficiently small, with bound depending on Cq: the fc < part of ()4.7.13p 
yields 

OiM < 0(0') + ^Te(0^-_0o) + 0(e'/'), 
so since 4> < and ag < 1 (for e small), 

0(0o)e(^,-) 

Contradiction! So (j)' cannot exist; 6{(f)) ^ (tt, ^ — agcj)) for any G [0, 0]. 

Now assume q{(t)o) 7^ and 0{(j)o) e (f , tt] for some 0o (see Figure [T7| . Let (pi g 
[0,(/)o) be maximal so that q{(j)i) = or 9{<j)i) ^ (f ,7r]. Again, p{<p = 0) = 0, so 
such a 01 must exist. For £ (0i, 0o], A;(0) > > p(0). 

The p < part of (|4.7.12p (in reverse direction) shows that q(0) > for all 
£ [01, 0o], so g(0i) = is not possible. 

For sufficiently small e, using (|4.8.6p . (|4.7.10p yields p^{4)) > for those 0. Thus 
p(0i) < p(0o), so 6'(0i) = f is not possible either. 

Only 0(01 ) = TT remains: but then p(0i) < = /e(0i), so (UTTTTT]) yields fc0(0i) > 
0. Therefore 0(0i — Si) ^ n + 62 for some small 61 £ (0, 0i) and 62 > 0. This is 
in the sector we have already ruled out — contradiction. The proof of (|4.8.3p 
is complete. 

Define 

q, p > 0, 



Q 



\k\, p<0. 
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Figure 17: The second derivative inequality (|4.7.13p . combined with the wall 
condition p = 0, rules out a sector of X4>i values. 

By (|4.8.3p . for p < necessarily k < pcoi{ae<j>) < 0, so 

= '^ =^^' -(yep < -(T6ifctan(tT6i0) = (tb ta.ii{ag (j))Q. 

For p > 0, (|4.7.12p yields Q^.^j < 0{Q'^). Altogether < 0{Q) (where O 

is with respect to IQj | 0), so integrating from <j> backwards yields 

Q(0)<O(Q(0)) (</'e[Oj]). 

Using \p\ < \k \ tan{ag(j>) for p < again, we get 

c-V-ig((/.) < c-V-iO(<7(0)) < Cptc'/^ (0 e [0,0]), 

for sufficiently large Cpt, and for e sufficiently small with bound depending on 
Cpt only. This imphes the strict form of (|4.4.22p . 

The strict form of (|4.4.23[) is immediate from = 1 — e, in the form 

xl^-\/{l-e)c^-x't 
(note that (|4.4.23p (nonstrict) fixes the sign). Here 6pn is sufficiently small. 



Density, velocity Finally, we obtain (|4.8.ip : (|4.4.6p combined with x = 
on the shock yields 

for some constant C. Integrating (|4.4.22|) along Pp we obtain 

niax|x-X^(G)l <CV/2 

Pr 
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for some other constant C . Combined with (|4.4.22p and = 1 — e this imphes 
(|4.8.1|) . for C, C", Cpp depending only on Cpt- 



(j4.8.2p is shown in the same manner. 



Left arc Pl can be discussed in the same fashion, by noticing that in the L 
picture (Figure [T^ it has the same properties as Pr, except for the wall not 
being horizontal which does not matter. Note that the mirror reflection in going 
back to the R picture reverses the direction of and therefore changes 9 to —9 
(see g;H3D). □ 



4.9 Corners moving along arcs 



Since our shock is a free boundary, we cannot be sure where the shock-arc corner 
is located. We study the behaviour of shock normal and downstream data when 
keeping the upstream data p„, Vu fixed, imposing Ld = \J\ — t. and restricting 
the shock location to be on Pr. (This is different from Proposition l2.6.2l where 
we imposed i"^ = instead of the location.) 

All the calculations in this section are done for the corner between S and Pr, 
but each has an analogous result for the left corner, using L coordinates (Figure 

[ED. 

We abbreviate w := 77^?. In this section refers to derivatives of upstream and 
downstream quantities as to is varied, while keeping ^r € Pr and maintaining 
the shock conditions and the parabolic boundary condition. 

We use dot notation = etc. on relations that hold only for u: = t]\. No d,^ may 
be taken of such relations. 

In addition we use the notation from Section t/„ — vi, z — Vx, v — ^Iil) etc. 
We required ^ e , so 

^=y^(l-e)4-7y2. 

Moreover 

n = (0,-l), f=(l,0) 

= -C, 4 = -ri, z^ = 7/, z^i = 
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Thus: 



d^n^ = -d^^l-in^)^ = d^n^ ^ 0, 

ny 

We use 

^ \z,\' = (1 - ^)c^ = (1 - ^)ci + ^^'^f''\ \zu\' ~ i.-.n 

7 + l + e(l-7) 7 + l + e(l-7) 



(PXTI) yields 



n.- ,2, , -2(7-1), 



9^.(1^-^1^) = ^dUlzul') + 0(e) = ^' ^ ' vl + 0{e). (4.9.2) 
7+1 7+1 

On the other hand, 

a.(|z,n = 2zld^{z-,) + 2z,^a.(z2) ^ -2Ca.(z,-) - 277a.(z,^), 

so 

^d^{zl) + 779^(z,^) = ^< + 0(6). (4.9.3) 
7 + 1 

(|2.3.2p can be restated 

= p^zl - p^z2 = Pu{zln- + zlny) - tt"! (^(p„) + (^^^ + 

Take d^j: 

=> = pu{n^d^z: + z^a^n^ + n^a^z^ + z^S^n^) - ^K^^dz-^l^) - dU\zdf))z^ 
- Pdirr'd^zl + z;^a^n^ + n^a^z^ + z^a^n^) 



2c^ 



^ .„(-ea.n^^ + 1) + ^g^^ - p.(-ea.n^ - d^zl) 



Pd c2(7 + l) 



^< + a.z,^ • ^ - (4.9.4) 



" ?7-i'm c2(7 + l) 
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using pu/pd = z^/z^ = 77/(77 - <). Finally, z*„ = z* yields 
so take du,: 

= - zl) + Tf{d^zl ~ d^zl) - d^ny{zl - zl) - ny{d^zl ~ d^z^) 

= -vld^n^ + d^z"^ - 7]/^ 

-vld^n'' + d^z^ = T^/i (4.9.5) 

(|4.9.3[) . (|4.9.4p and (|4.9.5p form a linear (nondegenerate) system for the three 
derivatives. Solution: 



h-iK.!-2(^i),-2,,,.,(,-„.i) ^ 

(7 + 1)(2?7 - -u^) 

Using i'^ = 2^ + '7: 



- ^<P^tS§^ + 0{e). (4.9.7) 



(7+l)(277-w2) 
Note that 

d^yy^ > (4.9.8) 
(for sufBciently small e) because ~ Vj < 0, rj — vy ^ z^ > 0. 
Transform (j4.9.3p to 

dMd) = ^jdM'd) + ^-T + Oie) (4.9.9) 

(no need to evaluate further; duj{n,^) is not needed). Finally: the counterclock- 
wise unit tangent for Pn is 

so 

P ^ X4' ^ i^Xt ^ rz ■ tp ^ ■ z ^ £,zyj^ - -qzl- 

Thus 

5F77((7+ !).?(- 277;) -2c? , ^^^^ 



(7+ l)(277 -z;il) ^ 
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Using T] = =■ crcu, v^^ = + -q = i] - z'^ = [a - i")cu, as well as 



c3 



= + ^{iz:r - (4)^) = (i + ~ a'))ci 

we obtain a more convenient formula: 

2 + (4.9.10) 

Since cr = z^jcu < z^jcu = ^2 for any admissible shock, we can argue that 

FTT-^ 7 — TT^ + '-'^^1 

K + a- (7 + lK 



> 



(7+ IK 

(<5,, + (7 - 1)^) • ^^f^ + 0(e) (4.9.11) 



for some iJp,, > 0. (Note: L" + cr is uniformly bounded because the set of possible 
shock locations is bounded.) Also, 

(4 K)^ = (c^). '^'^=^-' (1 - e)5.(|z.n ^ + 0{e) > 

7 + i 

(4.9.12) 

for e sufficiently small. 

^ ■ ^< = ^i^O . 0(, (4..13) 

(The O term is uniform in 7 j 1 because (— cr/Zcg)^^^ has a (7 — 1)^^/^ singu- 
larity which is cancelled by the 7 — 1 numerator.) 



Extreme corner locations Now we study the behaviour of p, k, q and 8 (as 
introduced in Section [4. 7p for rj ~ if^. We use a superscript ± to indicate quan- 
tities evaluated for 77 = 77^; a superscript * indicates Ty^j. We omit superscripts if 
the choice is unimportant (e.g. if the difference incurs an 0(e^/^) term which is 
dominated by something else). This is the case for p^^, k^^ and other derivatives; 
we may conveniently evaluate them at ?7|j (note Vr~Vr — 0(e^/^)). 

p* =0 and k* = 0(e) by (fUTE^ (with r = ^T~^cr and = cr). Then 

P+ ^pMi-V*B) + 0{e), (4.9.14) 

k+ = k* + fc„(,7+ - 7fj,) + 0(e) ™ kM+ - ,7],) + 0(e), (4.9.15) 
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where p^^, kui are the vahies at r] = rj'^. Therefore 
We estimate this: 



(4.9.16) 



V 

CdVu 



1 14.9. 13> 
1 14.9. 101 



2 I Tnf , 7-1 rn 



{LI + a)^/Cu 



')]e + 0{e) 



2 I T'nl , 7-1 T 



7 - 1 eg; - 
7+1 c2 



(is + 



0(6) 



9 r2 



7-1 

7 + 1 c« ' 7 + 1 " 



(is + ^) 



0(e) 



0{e) 



2 C2+2:2^(|£|2 



7 + 1 
2 



'u?~\z,?) , 7-1 2 /^ + i«(-+^^Ti«)A2' 

(7 — 



7 + 1 



c5 

7-1 



7 + 1 



{{Llf-^'')] + 



7-\ 
7+ l' 



(is + 

^ + iS('^ + ^i")\2 



(is + f^) 



0(e) 
l+0(e) 



2cS((iS)^-l) 



(7 + l)2c2(LS + ^)2 
The last factor is positive: (|2.4.18p yields 

7 - In 



(7 + l)^' + (7-l)(i"r +0(6) 



(4.9.17) 



a < 



-l: 



7 + 1 " 7+1' 



so 



2 - (7 + 1)<J' + (7 - l)(i^)' > ^^^^iK 1)'- 

7+1 

LJ]; — 1 is uniformly positive since the corner shocks allowed by (|4.4.6p are uni- 
formly not vanishing, for e sufhciently small. All other factors are trivially 
positive (note L" > l,a). The right-hand side of (|4.9.17p is positive, so 



and therefore 



j4.9.16t _ 

g+ < - 



V 

CRVJ 



ivi - Vr) + 0(e) 



(4.9.18) 



(4.9.19) 
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Extreme 6 value 

tanf- - e+) - ^ ™^ ^ + Ofe) ™^ /t+T + -^"(^ + ^-^") 
On the other hand: by convexity of tan on [0, 7r/2), with agcf) < (f) < 7r/2, 



tan(CTe(/') < ere tan0 = cre^ + 0(e) = \/^^ • 77— + 0(e) 

V 7 + 1 4/c« 

Then 

tan(ge0) ^ 7-1 + x 

tan(| - 0+) - 7 + 1 • _^ + + ^L^) 

The right hand side is < 1, for sufficiently small e, if and only if 

(7 - IML- + <j)<2 + L-{{j+ l)a + (7 - 

Since a = z^/cu < z^jcu = i," and > 1, this is always true. Therefore 

61+ e (0, 1 - (70 • 0) + 27rZ. (4.9.20) 



The result for follows from symmetry: 



€ (tt, Y - • + 2^2. (4.9.21) 



4.10 Corner bounds 

Proposition 4.10.1. For e sufficiently small: 

for any fixed point & T of IC, the lower bounds in ^J^J-^ strict: 

Proof. For 7 = 1: we may borrow (|4.8.8p which contradicts 

= ^ + PuiVR - V*r) + 0(e) = -Pu^e'/^ + 0{e) 
=0 

if e is sufficiently small, because > by (|4.9.10p . 

For 7 > 1: (|4.8.3|) contradicts (|4.9.21|) . for e > sufficiently small. 

Again the proof for the left corner is analogous: in L coordinates (Figure [T2|) 
Pl and ^i, ^2 have the same properties as Pr, S^r, except that £,bl is not 
on the horizontal axis which is irrelevant. □ 



74 




Figure 18: Left: {p, k) value curve (bold) from a k maximum on to the upper 
endpoint, for rjfi = iT^- The curve cannot go past 0+ because (|4.7.13p does not 
allow it to pass through the lower right quadrant "fast" enough. Right: {p, k) 
curve (bold) in opposite direction; the constraints limit the value of k maxima. 

To prove that ryj^ = r/^ is impossible, a more global argument is needed. 

Proposition 4.10.2. Consider rj^ — if^. There is ac so that for e sufficiently 
small, 

c(0)<c /or0e[O,0]. 

For 7 > 1, c satisfies 

c^^Jho + Cce+ {-(Tg/afY/'^q+ sin(6'+ + ag4>). (4.10.1) 
where C- is some constant, as in the O terms; for 7 = 1 loe may take any c > c. 

Proof. Consider 7 > 1. For A; < the result is trivial. Assume that there is 
a (/)o G [0, 4>) so that (j) ki'P) has a positive local maximum m. (j) = (see 
Figure ini)- For (/)o > we need ^^(^o) = which implies p((?!)o) = by (|4.7.11|) : 
for (/)o = we have p = anyway. By (|4.7.10p . fc > means > (if e is small 
enough, so that the O(fc^) = 0(e^/^fc) term is dominated by ugk). Therefore 
p > on ((/)o,(?f)2] for some 4>2 > 4>o- Take (j)2 € i<t'o,4'] maximal with this 
property. 

Assume that ip2 < <f>. Then p{4>2) = necessarily. 

(|4.7.12p for p > imphes q > on [(/)o, ^2], so is well-defined. Moreover k is 
strictly decreasing on [0o,(/)2], by (|4.7.1ip . 

Assume there is a i/ii G {(f)o,4'2] with k{(f)i) = 0, hence fc < on [i^i, 02]- (|4.7.13p 
(for fc < 0) integrated from 0i to e ^2] implies 

0{4>) > 0{(f>i) - + 0(e'/') > -|. 

=0 <1 <7r/2 
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In particular > 0{(j)2) > —tt/2 — contradiction to p{4>2) = 0. The assumption 
was wrong; necessarily fc > on [(/jq, 02]- 

(|4.7.13p (for fc > 0) implies that 9 is strictly decreasing on [(t)o,(f>2\ (for 0{q) = 
0(e^/^), i.e. e, sufficiently small). Then 6'(02) & (0, §) which contradicts p((/)2) = 
0. The assumption was wrong; therefore 4)2 = 0- 

By integrating (|4.7.13p (for fc > 0) from back to (/>: 

0{cl>)> 9($) + ae(4 -<!>) + 0{e'/^)^e+ + ae(4 -<!>) + 0{e'/^) (0G [0o,0]). 
Integrate ()4.7.12p (for p > 0) backwards: 

< q{l>) + 0(e) = q+ + 0(e) (</. S [(/-o J]). 

Then 



so 



= -f^ePi4>) = -crgq{(l))cosl 

> -agq+ cosi9+ + ae($ ~ 4)) + 0{e) (0 e [00,^]) 



fc((/)o) < fc((/)) + (709+ / cos (9+ + ae{<t> - 0))d<?i + 0(e) 



(j'*' sin 0+ — g+ ( sin 0+ — sin 



{9+ + ae{4-M))+0{e) 



< 



q+ sin {9+ + crg(j)) +0{e) 



Finally we obtain a bound for c: let (f>o be the global maximum point of k on 
[0,^], then 



supc(0)^ = ho + J^k{cf>) <ho + J^fc(0o) 
[0,0] V V 



<ho + ]J^1^ sin {9+ + ag(f>) + C-e 

for e > small enough. This is exactly the statement. □ 
Proposition 4.10.3. Again consider the case rjn = if^ . For 7 > 1 define 



- \ 2 

c 



a((/)) := Vl - e • I ( — ) - sin^ (/) -cos0), (4.10.2) 



1/2 



5(0) = \/l~7(c - Cfl) (4.10.3) 
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where c is as in Proposition (T^ ' for 7 = 1 take a = d = 0. For e sufficiently 
small, 

a ~ max 5,(0), 
4>e[o,4>] 

and ip + cannot have a local minimum (with respect to ^) on Pn U {£,bf/}- 



Proof. For 7 = 1: = 1 — e implies Xr > —\/T~~e ■ c = — r, so tpr > 0. By 
the Hopf lemma, this does not allow a local minimum of "0 at Pr. In ^br we 
argue that by Remark 14.4.51 we may consider the even reflection of across A 
which still satisfies the same equation, so the Hopf lemma also rules out a local 
minimum of ijj in ^br- 

For 7 > 1: 



' ' / /IT , _ / ^ ±^^ t |,pg ^ 



Call the right-hand side / and take of it: 



2(/ - cos(0)) sin - ^f^^^B cos (0+ + fTe(0 - 0)) 



2./ 



/ > by (|4.7.8p and r = — e • c^f, so the denominator is positive. By 
(|4.9.20|1 . cos(0+ +(To{4)— (p)) > 0, so the numerator is negative unless / > cos(f>. 
In that case = f + 0{e^/^), because / is 0{e^/^) due to q+ = ©(e^/^) and 

ha- c\ = 0(e) (by gZll))- But < < = arctan ^ < f - (5, with (5 > 
uniformly in e j 0, so for sufficiently small e there is a contradiction. Thus / is 
decreasing in 0, so it attains its maximum in </> = 0. 

For a local minimum (with respect to 17) of -0 + on Pr we need 

= (V' + 0^)0 = + 0^0 = V'0 - a'' sin 0, 

so "i/jt = r~^ipc/> = asin(/) there. In ^5^7 this still holds (minimum or not) because 
we have V'^ = = asin0. Another minimum condition is 



L — 1 — e 

> {-ilj + aS^)r ~ tpr + aS^r ~ r — J (l — e) — ijj'f + a COS I 



= \J\ — e ■ Cr — \j (?{\ — e) — sin^ + acos(/). 

This is equivalent (by squaring to eliminate the root and solving a quadratic 
inequality for a) to 
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But a > a{(j>) which, by Proposition 14.10."^ is greater than the right-hand side. 
Contradiction! □ 

Lemma 4.10.4. 

a< (4.10.4) 

for some da > (depending continuously on \), and e > sufficiently small. 



Proof. For 7 = 1 this is trivial since a = and v^j < < ^j?. 

For 7 > 1: (liTS)) (with = c\{l - e)) means ft-o = c|j(l + 0(e)). Taylor 
expand the square root in (|4.10.ip around h^: 

c = v/^+ 0(e) + -^q+i-ajaff'^ sin {6+ + ag'4>) + 0(e) 

Cfl + ^(-a<,/a/)i/2 sin {9+ + ag-^) + 0(e) 

il4 9.19[ _ J* 

■< c« + (^A - ^Jj)^(-'Tg/a/)i/2 sin (^?+ + ag-^) + 0(e) 



<CR+ (??;^ - r?|j) ( ^ - 



for some 5a > because {—(Jg/afY^'^ ~ 2^^^ + 0(e) < 2 for e > small 

enough and because sin (61+ + (70 0)) < 1 by (|4.9.20p . Use (|4.10.3p to get (|4.10.4p . 

□ 

Proposition 4.10.5. Consider a as in Ii4.10.3\ ). For e sufficiently small, the 
shock through with upstream data vj and pj and tangent (1, -^jw) "^d ^■ 



Proof. The shock through with tangent (1,0) is the R shock where = by 
construction, for the new slope and location differs from by (1) moving the 
shock up to from while keeping it horizontal, and (2) rotating it while 
holding it in to make its slope the above. Both of these changes are 0(e^/^), 
so it is sufficient to consider a first-order expansion using known derivatives. As 
before, we use uniformity in e J, 0. 

— u^, (T = ^ . 71 = —T])- '^d changes by 
V*R) + Oie). 



For (2), we use ([2XT|) with n = (0, -1): 

dfji'l = (dpvd) ■ n^n^ + {dfjVd) 
=-1 



For (1), we use p.4.19p (note ~ v ■ n = 
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Here (3 = arctan^^ = 0(e-'^/^). So changes by 



i l4.10.4l i / Bv-f^X fr, 

K ('-If (^^-""^-"a + w- 



(1) and (2) combined: the change is 

=ei/2 

By (|2.4.17p the first factor is > 2/(7 + 1). Therefore > for tjr = rj^, for e 
sufficiently small. □ 

Proposition 4.10.6. Let x & ^ be a fixed point o//C. For C^i sufficiently large 
and for e > sufficiently small, the upper part of J^.^.fip is strict: 

Vc<Vc {C = L,R). 

Proof. Let a be defined as in (|4.10.3|) . As shown in Proposition 14. 10."51 -0 + 0^ 
cannot have a local minimum at U {^br}- For rjji > rj^, we have (-0 + 0^)2 = 
V'2 > in by (|4.9.8p (for sufficiently small e), so the minimum cannot be 
in either (note that the domain locally contains the ray downward from the 
corner). 

On the shock: ^(j + = + a^, so 

J ^ ^ a wfsi + a 



(l + .s2)i/2 (l + s2)i/2- 
For a local minimum at the shock we need dt{ip + a^) = 0, so 



A global minimum, in particular < '4'{£,r)+cI'£,r, additionally requires that (as 
well as the rest of the shock) is on or below the tangent through the minimum 
point, because and thus + are decreasing in 77. Proposition 14.10.51 
shows that the shock through with that tangent has > for rjj^ ^ Vr- In 
the minimum point the tangent has same slope but is at least as high, so the 
shock speed is at least as high, so = -02 > is at least as high, in particular 
> too. But that contradicts a minimum (the ray vertically downwards from 
any shock point is locally contained in Q). Hence '0 + Q^C cannot have a global 
minimum on the shock. 
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ipn = on A contradicts a minimum on A (by the Hopf lemma) . 
The equation (|2.2.5p yields 

(a^ is linear), so the classical strong maximum principle rules out a minimum 
in the interior (unless -0 + is constant, which means we are looking at the 
unperturbed solution which has Vb. = Vr < Vr)- 

On Pl: for 7 = 1 we argue that ipr < on (in fact ipr < except in 
the unperturbed case rjl = rj'^), as in in the first paragraph of the proof of 
Proposition 14. 10.^ so the Hopf lemma rules out a maximum at Pl U {£,bl}- In 

either ijl < vl < V*r < Vr so that ^(6.) > "0(6?.) = V'(CA) due to Tp = 
on S and V'| = > i/',^, or rjL > r]*^ so that V?? > in ^l, by the same analysis 
as for ^ij, ruling out a minimum. 

For 7 > 1: 

+ «0« = ^« + 0(eV2) + 0(,i/2). 

For C,, sufficiently large and e sufficiently small^the right-hand side is negative: 
u£ = for 772 = Vr, it is strictly decreasing in (Proposition 12.6.^ . and Crj 
bounds away from rj'^ (see (|4.3.4I) ). Again, no minimum is possible. 

Conclusion: no ip minimum anywhere — contradiction! 

The argument for rj^ is similar, using the L picture fFigure Il2p . The wall 
still passes through the origin, so the wall boundary condition Xn = implies 
tpn = 0. For 7=1 this implies that ip cannot have minima at the wall (Hopf 
lemma). For 7 > 1 wc have to modify the argument at the wall since it is no 
longer horizontal. A minimum of ip + aS, requires [tp + a^)„ < (where n points 
out of il), so a^ri < 0, so ^„ < (by a > 0). But the wall slope is positive, so 
> (for outward n) — contradiction. 

The remaining arguments are as in R coordinates. □ 



4.11 Density bounds and shock strength 

Proposition 4.11.1. If e and dp are sufficiently small (bounds depending only 
on Cpt), then for any fixed point ip & of K,, the inequality l4-i-9^ is strict. 

Proof. By Proposition I4.5.2[ ?/; and hence s are analytic. Thus we may use 
Proposition [3?2Tl] which rules out minima of p in 57 and (using Remark 14. 4. Sp at 
A, unless ip and s are linear so that p = pR is constant and (|4.4.9p is trivially 
strict. However, p G C(ri), so it must attain a global minimum somewhere. 

^^We need Crj > 0; otherwise the upper bound on e would become zero as r)£ | r?|j. 
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p minimum 




Figure 19: In a p = p_ minimum at S, the shock tangent must be above ^ and 
^/j, so by (|1.3.1[) the minimum is separated from Bcj{vi). Therefore p- ^ pi- 

Case 1: p attains its global minimmn on S (excluding ^l,r)- Then by Propo- 
sition [XTHl the rest of the shock, including the corners must be below 
the tangent in that point (see Figure [TO]) . We assumed in Definition 14.3.11 that 
the line through ^2 and does not touch or intersect the circle with radius c/ 
centered in vj. Hence if e is sufficiently small (depending on rj * ), then by (|4.4.6p 
the shock tangent must have positive distance from the circle as well. Thus the 
global minimum of p must be bounded below away from pi. 

Case 2: p attains its global minimum on Pl or Pr . By (j4.8.ip we know p up 
to 0(ei/2). 

Combining both (nonexclusive) cases, we see that for sufficiently small Sp and 
e, depending continuously on Cpt, A, if^, ()4.4.9p is strict. □ 

4.12 Velocity and shock normal control 

Proposition 4.12.1. If Ssa is sufficiently small, if Cyx,C'vL cire sufficiently 
large (bounds depending only on Cpt), if Csn is sufficiently large (bound depend- 
ing only on Cyx,CyL), if S^y and e are sufficiently small (Syy bound depending 
only on Sp, Csn', e bound depending only on Csn), ond if 5cc is sufficiently small, 
then for any fixed point if^ £ J- of K, the inequalities ^.4-4\ ), H-4-^4V y i4-4-25\ ), 
i4.4.26\) , 14.4.27^ , <^re strict. Moreover 

\Xt\>S^t (onSnBs,iic),C = L,R) (4.12.1) 

for some constants (5^t, Sd > 0- 

Proof. 1. For horizontal velocity (j4.4.24p : assume that attains a positive 
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global maximum (with respect to Q) in a point on S. Since vj = (0, v^) 

with Vj < 0, this means < (because < 0), i.e. si(^o) < 0. 

Si(^o) can be expressed as a continuous function of v^{^q) and ^o- The set 
of possible is contained in the set of possible shock locations which is 
pre-compact. Therefore if has a maximum — C^j.e^/^ in some ^ '5, 
then 

si(Co) < (4.12.2) 

where Csi = Csi{Cyx) > is uniformly increasing in C^x • 

Since S and -0 are analytic (Proposition 14. 5T2|) . we can apply Proposition 
13.5.11 with w — (1, 0). For a constant-state solution (|4.4.24p is immediate. 

Otherwise p.S.ip and (|3.5.2p are satisfied, < means < 0, so by 

([5X^ K < 0, i.e. sii > 0. 

Now we can use a similar argument as for the density in Proposition 
14.11.11 and Proposition 13.7.11 (see Figures [TUl and [TO]) : sn > imphes 
< si(5o) for £, < £,0 near ^o- On the other hand, for £l = CI "we have 
si(Cl) > by construction of the L shock in Section [Ol so for arbitrary 
5l (satisfying (|4.4.6p ) the continuous dependence of si(Cl) on rjL and by 
(14X6)1 shows si{£_l) > -C2e^^^ for some constant C2 independent of e. If 
we pick Cvx so large that Csi > C2, then si(Cl) > si(Co) by (|4.12.2p for 
any possible location of the left corner. 

Therefore we can pick G {£,lAo) niaximal so that Si(Ca) = Si(Co)- Then 
SiiO < si{£.o) for C e (Ca,Co), so by integration 

But that means the shock tangent in is parallel to the one in but 
higher^ so cr := ^ • n is smaller in By (|2.4.19p . that means is smaller 
in whereas u* is the same (parallel tangents), < 0, so is bigger 
in Contradiction — we assumed that we have a global maximum of 
in Co- 

Propositions 13.3.11 and 13.4.11 rule out local maxima of in fl and on A, 
where we use that x is analytic and that w is not vertical, i.e. not normal 
to the wall. 

On Pl U Pr we can use (|4.8.2p . increasing Cvx to > Cpv if necessary (this 
makes Cvx depend on Cpt as well). Now (|4.4.24p is strict. 

2. In any point (z S, si{C) is a function of v^, with sgnsi = — sgnii^. 
si(C) is continuous in ^ and v^, and the set of possible shock locations ^ 
is pre-compact, so (|4.4.24p implies 

sup/(r?, nfl) < Cs„e^^^ 

where Csn = CsniCvx)- 
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3. The arguments for (|4.4.26p are anafogous to those for (|4.4.24p : the trans- 
formatfon from R to L coordinates (Figure [12]) turns n-j^ into nj^ = (1,0) 
(and vice versa). The wall is never vertical in L coordinates, so Proposition 
13.4.11 still applies to v^. The other arguments arc as before. 

Moreover (|4.4.26[) implies 

supZ(nL,n) < C5„e^/^, 
where Csn = Cs7i{Cvx,Cvl) now. (|4.4.27p is strict with these choices. 

4. 1)2' = on A; = 0{e^/^) on Pl U Pr by (|4.8.2p . has no cxtrema in fl 
f Proposition 13.5.1]) : to show (|4.4.25p it remains to discuss S. 

()4.4.9p is in particular a lower bound for the shock strength, so \v2 — f"| 
is bounded away from 0. (|4.4.27p bounds the shock normals away from 
horizontal. Both combined imply ()4.4.25p is strict if S^y > and e > are 
chosen small enough, with upper bound on Syy depending on Sp and Csn, 
and upper bound on e depending on Csn only. 

5. The shock normal bounds also imply (|4.4.7p is strict, for Scc > and 
e > sufficiently small(er), with e bound depending only on Csn- Here 
we use (|4.3.2p : (|4.4.27p shows shock tangents are between and n-^ (up 
to 0(e^/^), and these are bounded away from arc tangents. 

6. Since the left corner is above the wall, the shock normal boimds imply 
(|4.4.4p is strict, for e sufficiently small(er), with bound depending only on 
Csm and for 5s a sufficiently small. 

7. Near each corner the shock normal bound bounds n away from the ^ 
direction, so \xi\ > S^t and therefore ()4.12.ip for some S^t- 

□ 

4.13 Fixed points 

Proposition 4.13.1. For So sufficiently small, with hounds depending only on 
Sp and Cl, for Cd resp. Sd sufficiently large resp. small, with bounds depending 
only on Sp and Cl , and for e sufficiently small, with bounds depending only on 
Cpt, Cl and Sp: 

If X ^ ^ is a fixed point of IC, then and {^.^.29^ are strict. 

Proof. First we check (|4.4.28p . 

1. For the wall boundary operator (|4.4.19p . the strict inequality (|4.4.28p is 
obvious. 
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2. The parabolic boundary operator (|4.4.17p has p derivative = Vx. We 
use (|4.4.22p with sufficiently small e, depending on Cpt, to obtain (|4.4.28p 
strictly (with 5o — j for example). 



3. For the shock boundary operator (|4.4.18p 
vi — Vib , ^ „ - 



9} - Pd - c-Vx^)^L^ - \v, - V^l- (l - (^^)^) • (PVX - P.V/) 

\l)T — \/lh\ ^ I?l7- — ' 

(4.13.1) 



We exploit that for a fixed point we have tf; = ijj = ij/ at the shock, so 

n ~ — — J 

\vi - VV'I |Vx' - Vxl 

where n is the downstream shock normal. Note that 

|Vx^ - Vxl = xi-xn 

because xl = Xt and xi> Xn- 

g}; = p(l - c-2Vx')?i - \vi - V7/.|-irr- (pVx - P/Vx^) 
p{l - c"2Vx^)n- (xfi - XrO^^Cp- P/)Xt** 
'^p((l-c-^Vx^)n-(x„)-^Xtt) 

= /5((1 - C-2;^2)^ _ Xt{{Xn)-^ + c'^XnY) (4.13.2) 



=0 

= p{l - c-'xl) > 5'o 

for some So' depending only on Sp, because (|4.4.9p bounds downstream p 
away from pj, and that means 1 — c^^Xn ~ 1 ~ i-^d)'^ lower-bounded 
away from 0, by the shock analysis in Section [531 1 5^1 is easily bounded 
from above, using (|4.4.2ip . so (|4.4.28p is strict for a sufficiently small So, 
depending only on Sp and Cl- 

Now we check (|4.4.29p . For the parabolic- wall corners it is trivial: on Pl or 
Pri 9p = Vx is almost normal, by (|4.4.22p using a sufficiently small e (bound 
depending only on Cpt); on A the vector g^^ is normal; the corners enclose an 
angle exactly 7r/2, so the g^ directions are independent (almost orthogonal). 
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For shock-parabolic corners: on Pl resp. Pr, again gt = V^. On S, use (|4.13.2p . 
We normalize both derivative vectors and consider their cross product: 



^/ii-c-'xir+xUix'n)-'+c--'xny 

^ Xt(l +Xn/X») 



(4.13.3) 



Denominator: |Vx| > \/l — e ■ c > c/2 for e < ^. The square-root is lower- 
bounded by 1 — c^^Xn > > for some S depending only on Sp (see above). 
(|4.12.ip bounds the numerator away from 0. 

Hence ()5.1.9p is strict, for d sufficiently large, Sd and then e sufficiently small, 
depending only and continuously on Sp. □ 

Proposition 4.13.2. If the constants in Iji4.3.3\ ) in Definition \4.4-3\ are chosen 
sufficiently small resp. large: 

for any A G A, IC\ cannot have fixed points on T\ — T\. 

Proof. Let x € be a fixed point of /C. We show that every inequality in the 
definition of T is strict, so x G 

(|4.4.2ip and (|4.4.2p are strict by Proposition 14.5.21 

(|4.4.4p is strict by Proposition l4.12.1l 

(|4.4.9p is strict by Proposition 14. 11. II 

A fixed point satisfies t/' = ^, so Wi/; — ^\\ — rj{ip) > cannot be true. (|4.4.12p 
is strict. 

(|4.4.10p strict is provided by Proposition 14.6.11 

Due to Proposition 14. 6.11 = 1 — e on each point of Pl U Pa, so we are in the 
situation of Section |47] etc. Proposition l4.8.1l shows that (|4.4.22p and (|4.4.23p 
are strict. 

(|4.4.24p is strict by Proposition l4.12.1l 
(|4.4.25p is strict by Proposition l4.12.1l 
(|4.4.26p is strict by Proposition l4.12.1l 

Propositions 14.10.1] and 14. lOTBl rule out t]l = ril± S~^e and rm = ri*j^ ± S~^e if S 
is small enough, so (|4.4.6p is strict. 

(|4.4.7p is strict by Proposition l4.12.1l 
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(|4.4.27p is strict by Proposition l4.12.1l 

Proposition 14. 13. II shows that (|4.4.28p and (|4.4.29p are strict. 
All inequalities are strict. 



□ 



4.14 Leray-Schauder degree 

We determine the Leray-Schauder degree of /C on for a particular choice of 
parameters A: a straight horizontal shock (yy^ = r]^), with 7=1 (see Figure 
[^U]) . This problem is simple enough to compute the degree precisely, although 
the discussion is still difficult, especially due to the free boundary. 

Another option is to introduce a homotopy to an even simpler problem (linear, 
fixed boundary, ...). But it would be necessary to prove a new set of a priori 
estimates for a family of arbitrary, unphysical problems, hence a lot of work 
without useful sideeffects. 

Proposition 4.14.1. For sufficiently small e: 

For 7=1 and if^ = 7]^, there are no fixed points of IC in J- other than the 
unperturbed solution. 

Remark 4.14.2. This result makes essential use of 7 = 1 as well as the con- 
siderable simplifications from using the constant-state solution. At the time of 
writing we do not know a way of proving uniqueness for 7 > 1 or Mj < 00; this 
is why we choose to argue by homotopy to 7 = 1, jy^ = 77^. 

Proof of Proposition \4.14-1\ Let ip be the unperturbed solution and s the cor- 
responding shock (s(^) = 77^ for all ^). Consider another fixed point ip for the 
same case 7 = 1 and ry^ = rj'^- Let s be its shock. 

Consider ^ coordinates as defined by -0 = ip{a, C,) (see Definition 14. 4. 3p . Then 
= (c^/ - Vx^) : V^V = [c^I - Vx^) : V^(?/' - V') (in ^) 

because is constant. The classical wcalP°l maximum principle |GT83| Theorem 
10.1] implies that ~ must attain its global minimum and maximum on the 
boundary. This excludes the bottom boundary because we can make it interior 
by reflection (Remark 14. 4. 5p . 

Assume ip has an extremum (with respect to fJ) in ^ G U Pr. Then ipt = 
there, so Xt = because ^ ■ t ~ on Pr and Pl for r/l — ifj^. So b\F^ 

^"The strong version applies, but is not needed 
^^This argument breaks down for 7 > 1. 
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Figure 20: The unperturbed solution: a straight shock. This is the asymptotic 
limit for M/ 1 oo while holding Mj sinr constant. In this limit, the wall corner 
moves to ^ = (— oo,0). 



|VxP = (1 - e)c^ we get Xn = -(1 - eY^^c (the sign is fixed by ()4.4.23p ) and 
finally 

'4'n ^ Xn + S. ■ n = -cVl - e + cVl - e = 0. 

Again there is a contradiction to the Hopf lemma. A reflection argument (Re- 
mark H74.5p also rules out extrema (with respect to Q) in S,bl and S,br- 

Assume ip has a maximum > ip on S. By ip = on S, strictly decreasing 
in ?7, this corresponds to a minimum of s smaller than s = ifj^. Then 8^=0 
in that point, so graphs is horizontal there, like graphs. Since it is lower than 
graphs the shock S is weaker there, so by (|2.4.19p we have ^2 < 0. But this is 
incompatible with a maximum. 

Analogously a minimum < ip is ruled out. 

Assume 1/' ha-s a maximum > 1/; in ^r. Then ip ~ on S, decreasing in 
77 means s(^_r,) < Tyjj, so ^2(Cfl) < by (|4.9.8p . for sufhciently small e. This 
contradicts a maximum of ip, because the negative vertical direction from is 
contained in the domain. After reversing some signs and inequalities we rule 
out a minimum < ip in the same ways. Analogous arguments apply to ^l- 

Altogether we must have tp = ip. □ 



Proposition 4.14.3. The unperturbed solution with 7=1 has index l ^ as 
a fixed-point of K, from Definition \4^.4.9\ 

Proof. We use |Zei86| Proposition 14.5]. Since JC is compact fProposition l4.5.3|) . 
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we have to show that / — dJC/dip{tl}) has trivial kernel, where dlC/dip{il)) is the 
Frechet derivative K. a.t ip = ijj. If this is true, then the index of ■)/) as a fixed 
point of /C is ±1. 

We consider ?/' = = V' in (? = (cr, () £ [0, 1] coordinates (from Definition [4A3|). 
We consider first variations ip' of ip in these coordinates: i.e. consider ip + tip' 
for t e IR and evaluate dt a.t t = 0. Let /C('0)' etc. denote the resulting first 
variations of other objects. 

Assume that ip' ~ dlC/dip{^)ip' (which is /C('i/;)')- We have to show tp' = 0. 

tp' = IC{ip') implies that the variations of the n- ^ transforms defined by ip 
and JC{ip) (see Definition 14. 4. 3p are identical as well. We may write (^)' without 
distinction. Moreover ip' = ip' . 

The following relations are all meant to hold for t/i ~ V ~ V' and ip' = ip' only. 

We emphasize here that the variation ip' is taken in a coordinates. The vari- 
ation of X = "0 — is ''T-ot (necessarily) ip' — because varies as well. 
We deliberately vary tp, not x> because VV' = in all coordinates, allowing 
considerable simplifications. 

In particular: 

(V^-^)' = (VpVffV^)' 




=0 



= (E|J)'vp-+i:|^(vp)' 



=0 

+ (V|(7)'v|VV^-t7 + Vp(V|V;)'V^-^ + WfaVitPiV^-a)' 

=0 =0 

The same relations hold with ip instead of tp. 
Frechet derivative of the interior equation: 

= (c^/ - (V^-x)^)' : V|V3 + {c^I (V^-x)2) : (V|^)' 

=0 
2„/,' 
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The resulting right-hand side is a linear elliptic operator without zeroth-order 
term, applied to ip' . The classical maximum principle shows that ip' cannot 
have a minimum in the interior. 



On the paraboHc arcs, 
linearizes to 



|V4I 



= Vrx • (v.-x)' = (v.-v- - • ((v.-V-)' - iO') 



v.; = 



= (r2)'=0 



(4.14.1) 



Here we use that variation of -0 may move ^ but keep it on Pr. (|4.14.ip does 
not admit any extrema of "0', by the Hopf lemma. 

By reflection across A we can also rule out extrema on A, by applying the 
arguments for Pl and Pr resp. Vl there. 



Shock: 0; = = V'^(0, 0) + wfr/, so 



on shock 



(4.14.2) 



Moreover 



( ^7 - Vg-V' ^, _ -1 



I 



\ 



\ 



Wi ~ V/01 



(4.14.3) 
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Take the Frechet derivative of (|4.4.18p : 



o=((pV^^x-P/V^-/)-— — 



= pc-\- X' - V-x • (Vrx)')Xn + P(V.-V' - 0' • « + ■ n + (pV-x - P/V^) ' (- -^t)' 

^ ^ ^ \vi — y ^\ 

=0 

(pV,-x - P/Vrx^) • 

= -pc-^Z + P(l - c-\l){^\ - pc-\nXM')t + (p - pi)v' - jJ-^Xt^t 

\vi - V|^| 

= p(i - c-\i){i,% + (^^s^ - pc"'x«)xt(^')t + (p - piW - pc~'m 

\\vi~V^\ I 
|4.14.2[ _2 , ( Pi-P -2 \ f I,-. , (P- PI -2^ ,/ 

= P(l-C X„)(V')n+ (^r^— T-PC XnjXHV' )t + (— J^y PC )^ 



^° " ^0 ' ^ <0 ^ 

(4.14.4) 



Now consider the case of a positive maximum of ip on the shock. Then {ip')t = 
{tp')i = which also implies (n)' = by (|4.14.3p . (This is natural because then 
we perturb the shock in a way that keeps it horizontal in that point, so the 
normal docs not change in first order.) The two remaining terms: 

= -P(l - c-\jm2 + - pc-')^P' 

> „ ' w| 



<0 



<0 



For a positive local maximum of ip' we need {'ip')2 > 0, but this is incompatible 
with the previous equation. In the same way a negative local minimum is ruled 
out. 

Finally consider a positive maximum of ip' in The parabolic condition 
(j4.14.ip means that V|V'' must be tangential to Pr in the corner; for a maximum 
we need that it points counterclockwise, hence into the upper left quadrant, 
unless it is zero. On the other hand the vector consisting of the {ip')2 and 
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coefficients in (|4.14.4p points into the lower right quadrant (because xi < in 
the right corner), so the scalar product is < 0. But the coefficient of is < 0, 
so the right-hand side of (|4.14.4p is < — contradiction! 

The same argument, with certain signs and inequalities reversed, rules out a 
positive maximum in the left corner as well as a negative minimum in either 
corner. 

We have ruled out that ij)' has a positive maximum or negative minimum any- 
where in the domain. Therefore "tp' = 0, which is precisely what we had to 
show. □ 

4.15 Existence of fixed points 

Proposition 4.15.1. For sufficiently small resp. large constants in /C 
has a fixed point for all X G A. 

Proof. Let t G [0, 1] X{t) be any path in A. We have shown in Proposition 
14X81 that U := Ute[o,i] (W ^ ^Mt)) is open in [0,1] x C|'"([0, 1]^), that its 
closure is contained in W := Ute[o i] ({*} ^ ^\{t))^ so 

dU CV -W-U ^ y (Wx (:fA(t)-^A(t)))• 
te[04] 

By Proposition 14.5.31 /C is a continuous and compact map on W . Finally, by 
Proposition 14. 13."2l we know JC cannot have fixed points on V , hence on dU . U 
is bounded. Thus we may apply the property (D4*) in [Zei86[ Section 13.6] to 
argue that the Leray-Schauder index of A^A(t) on •^A(t) is constant in t e [0, 1]. 

For A(0) = (7, '72) = (l:*?]?); Proposition 14. 14.31 shows that the Leray-Schauder 
degree oi JCx is i 7^ 0. Since A is path-connected (Lemma 14. 3. 3p . there is a path 
to any other A = A(l). So K,\ has degree i 7^ for all A £ A. Nonzero degree 
requires at leasll^ one fixed point. □ 

4.16 Construction of the entire flow 

Proof of Theorem[l[ For all pj, ci G (0, 00), AIj G (—00, 0), for each choice of 7, 
77* and for all sufficiently small e, a separate A is defined. For sufficiently small 
constants in (|4.3.3p . Proposition 14. 1 5 .Tl vields fixed points V' for all A G A. Note 
that there is no lower bound on e, except that a, (3 and some other constants 

^■^We expect that there is only one, but we do not need this fact. For all we prove, there 
could be more, e.g. two with index t and one with index — t (which would still yield total L-S 
degree t). 
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Figure 21: Extension from to entire domain 

may deteriorate as e j 0. By Definition I4.4.3( Remark 14.4.41 Proposition 14.8.11 
and ()4.5.3|1 . the fixed points satisfy 





= 


in 17(<^), 


(4.16.1) 






and 


(4.16.2) 


\P- Pel 




and 


(4.16.3) 


VV' - vc\ 


= 0(ei/2) 


onF^"' {C = L,R), 


(4.16.4) 


X 




and 


(4.16.5) 




= 


on S, 


(4.16.6) 


Xn 


= 


on A, 


(4.16.7) 




= 0(ei/2) 


(C-L,i?). 


(4.16.8) 



where the O constants are independent of e. For regularity we have 

ll^llco.i(fT<") ^^1' 

where d ~ d(X, p|'^ U P^^) > 0. 
for constants Ci and C2((i) independent of e. 

We extend each ip to a function ^/'('^^ defined on CW as shown in Figure [21) we 
use pL,VL in the region enclosed by wall, L shock and , we use pb,vb. in 
the region enclosed by R shock and , and pi,vi elsewhere. In each of the 
four regions, 4'^'^^ is a strong solution of self-similar potential flow, so we may 
multiply the divergence form (|2.2.3p with any test function 'd G C^(CW^) and 
integrate over each of the four regions to obtain a sum of boundary integrals. 
Set 

Dc := Pc^^^ - Pi^\ Ec := P^'^ - P*c^'\ {C = L, R). 



(4.16.9) 
(4.16.10) 
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These curves have length 0(e^/^), so the integrals over them can be neglected 
because ■0^'^) has e-uniformly bounded C^'^ norm in each region. The two inte- 
grals over S''^') cancel because of (|4.16.6p . The two integrals over Pj^^^ n P^^'"'^^ 
are 0(e^/^) due to (|4.16.3p and (|4.16.4p . and same for Pr. The integrals over 
each piece of the wall vanish because Xn = in all cases. Therefore: 

f p('^)Vx^'^ • Vi? - 2p(')i? d^"= 0(6^/2). (4.16.11) 
Jew 

By (|4.16.10p with a diagonalization argument, for every compact K C £iW — 

P]^^'' — Pfi^^ we can find a sequence (sk) i so that 0^'^'=) converges to -0'"^ in 
C"'^(A'). Moreover, the convergence is bounded in C°'^(CVF), so we can take 
e i in (|4.16.1ip to obtain 

f p(")Vx'°^ • Vt? - 2p('')z9 d|'= 0. (4.16.12) 
Jew 

In addition, (|4.16.5p . (|4.16.2p and (|4.16.9p show that 

e C{ZW) (4.16.13) 

Finally, by construction of ip'"'^\ 

p^°\sO,v^°\s^) pi,vi (s^oo) (4.16.14) 

for any ^ G £tW — {0} (note that any ray from the origin is either in W or enters 
and stays in the / region). 

(|4.16.12p . (|4.16.13p and (|4.16.14p show that 0(t, x) := V^"^^"^^) defines a solu- 
tion of (HXH), (fTXa . (fTXD . 

Taking 7 "f oo and ^7^(7) i ^l^il) obtain a solution for every 7 G [1,cxd) and 
Vl G {vI{i),V*r]- Note that 77^ T V*r for e i 0, by Definition gXH 

It remains to show that we have covered every tip shock that satisfies the condi- 
tions of Theorem [TJ Consider standard coordinates (Figure [12). By Proposition 
I2.6.2[ the shocks with some upstream velocity vi = {0,v^), Vj < 0, and density 
Pi and downstream velocity parallel to the wall are uniquely determined by rj^ . 
The choices with rj^ > rj'^ cannot intersect the horizontal axis left of ^2, as 
their shock slope is negative by Proposition 12.6.21 The choices with 77^ < 
are possible, but they are supersonic-subsonic at the intersection point (note 
that their other = 1 point must be above the axis, as their slope is positive 
and dBc^{vd,) is centered on the horizontal axis). Therefore, all supersonic- 
supersonic tip shocks have ry^ G (0, 77^] after passing from original to standard 
coordinates; there is only one such shock for each r;^ and we have constructed 
a solution for each that satisfies (|1.3.ip . □ 
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Remark 4.16.1. It remains to discuss whether there are tip shocks that satisfy 
the conditions of Theorem [TJ Proposition I2.6.2[ as discussed in the last para- 
graph of the proof of Thcorcm[Tl ahcady settles that question. By Lemma l4.3.2l 
for any Mf < 0, there is an interval ('7£(7), '7|j] of rj^ so that (|1.3.ip is satis- 
fied, ryf < 77|j, so there are some nontrivial solutions. For small \Mj\ we have 
rji^ > 0, so there are cases where (jl.3.ip is violated. But if \Mj\ is sufficiently 
large, then rjf^ = 0. In particular if Mj < —1, then Bc,{vi) is below the wall 
(in R coordinates), so (jl.S.ip is always satisfied. 

Remark 4.16.2. In addition to mere existence we obtain some structural infor- 
mation in the proof: 

1. The solution has the structure shown in FigureHl with L > 1 in the /, L, R 
regions, L < 1 in the elliptic region. 

2. The solution has constant density and velocity in each of the /, L, R re- 
gions. 

3. The solution is analytic everywhere except perhaps at the L, and fi, R 
interfaces. 

4. The shock is Lipschitz; it is straight between the /, L and /, R regions; it 
is analytic between / and the elliptic region (away from the corners). 

5. Density and velocity arc bounded. 

It is expected that density and velocity are at least continuous. However, the 
methods developed in this article yield boundedness everywhere, but continuity 
only away from U P^. Note that P^ , P^ can not be classical shocks with 
smooth data on each side, because the one-sided limit of L on the hyperbolic 
side of the parabolic arcs is = 1 everywhere (> 1 is needed for positive shock 
strength) . 

In some points, continuity can be obtained a posteriori. For example < 0, as 
implied by (|4.4.24p . yields continuity in £,br- 

Some other structural information: 

1. The possible (downstream) normals of the curved shock portion are be- 
tween the R and L shock normals (counterclockwise). 

2. The shock is admissible and does not vanish anywhere. 

3. Therefore, the shock is above the line connecting the arc-shock corners. 

4. In the elliptic region, the velocity tangential to the wall is between vf^ 
to v^. 

5. In the elliptic region, the density p is greater than pj. 

Additional information can be obtained from the inequalities in Definition l4.4.3l 
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Figure 22: Corner 

5 Appendix 

5.1 Regularity in 2D corners 

We obtain C^'" bounds in a corner neighbourhood from merely (with a uni- 
form C°'^ bound). Although our argument follows |Lie88a| Theorem 2.1] without 
major modification, we prefer to repeat it explicitly, because we need a uniformly 
large neighbourhood. 

We keep the proof simple by restricting ourselves to smooth coefficient functions 
and a priori rather than C*^'^ regularity. 

Proposition 5.1.1. 1. Consider a point xq G IR^ and polar coordinates {r, (j)) 
with origin in that point. Let R > 0. Consider two functions (p^^cj)^ : 
[0,R] [-f,f] so that 01 < (j)'^. For j = 1,2 let H be the curve 
{(r, 0J (r)) : r S (0,i?)}. Assume that H are C^'" curves (a e (0, 1] 
including the endpoints: 

< Cr. (5.1.1) 

Let I = {{R,(f>) : (j) € {(j)^ {R) , (j)"^ (R))} . Let U = {(r,0) : r G (0,i?), S 
((/>i(i?), 0^(i?))}. Assume there is a < ir so that 

sup cf>\r)<l, inf 0i(r)>-^. (5.1.2) 

Consider a function u £ C^{U) D C'^{U U U T^) n C^(U) with 

Ikllci(F) < Cu. (5.1.3) 

Assume that u satisfies 

a^'{x,u{x),Du{x))uij{x) =0 Va; e [/ (5.1.4) 
fwe use Einstein convention in this section) as well as boundary conditions 
g'^ix, u{x),Du{x)) =0 Va: e ^^ /c e {1, 2}. (5.1.5) 
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Here a'^ and g'^ , as functions of x,u,^u, are assume^^ to satisfy 

ll'^'^llci(I^xlRxlR2)i ll5''llci(t^xlRxlR2) < ^c- (5.1.6) 

Assume uniform ellipticity: there is a Se > 0, independent of x Cz U , 
z p e IR^ and ye IR^, so that 

a*^(f,z,p)yy" ><5e|yp. (5.1.7) 

Moreover, assume uniform non-degenerate obliqueness: there is a So > 0, 
independent of x £ F*^, 2 e IR and p £ IR^ so that 

g';{x,z,p)-n{x)>5o, (5.1.8) 

where n is the outer unit normal to F'^ (k = 1,2). 

Finally let Cd he independent of x £U , z £ IR and p £ IR^, so that 

\\G-^\\<Cd, (5.1.9) 

where G is the 2x2 matrix with kth column gp{x, z,p). 
Then there are R' £ {0, R) , X > and Cr < 00 so that 

\Vu{x)-Vu{xo)\<Cr\x-xo\^ forx£U. (5.1.10) 

R',X,Cr depend only on Cd, Se, Cc, 5o, Cu, Cr and R, except R' may 
also depend on a . 

2. Assume in addition that 

ll'-^llci'<'(FxiRxiR2) < C'r. (5.1.11) 
Then there are k > 0, Cr < 00, R' > so that 

ll-llc-(B.,(.„)nT7.{.o})<a. (5.1.12) 

A is as before. R' may have changed, but has the same dependencies. Cr 
may have changed and may depend on a now. k depends continuously and 
only on Cd, Se, Cc, So, C„, Cr, a, R. 

Proof. In this proof let C, Ci, S and Si (for numbers i) represent constants that 
may depend continuously and only on Cd, Se, Cc, So, C„, Cr, R (but not a). 
C, 6 may change from occurence to occurence. (5 > is meant to be small, 
|C| < 00 large. 

Constructing w: interior 



'Note that is defined for all arguments. See |Lie88a] for more general circumstances. 
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Define 

v{x) := g^{x, u{x), Du{x)) + AIw{x). 

Set 



u;(r,^) :=rV(0) /(0) 1 - 



with A e (0, ^] and S > to be determined. We take ji so small that f > \- 
In a given point of f/, rotate coordinates to be angular and readial and let 
^rr ^ ^rc/>^ ^(fxp -j-j^g corresponding components of a^K 

= {a''''\{\ - l)f + 2a''*(A - 1)/' + a'^'^(/" + A/))r^-2 

= ( a''''A(A - 1)/ + 2a'''^(A - 1) ■ fiBe-^'^ + a^^{X f - MB^e-^"^)) r^-^ 

< (e--«VS(2|a'-"*| - Ba'f"f') + Xa't"t')r^-^ 
We picl0 

B = (1 + 2sup la''"^!) / inf a-^-^ > 

so that 

a'^ Wy < ( - fiBe-^''' + Xa''"t')r^-^ 
For any ^ > we can pick a A = A(/i) > so small that 



Boundary condition for w 

Let i^, T be Cartesian coordinates that are normal resp. tangential in a given 
point on (then = 9„ and dr = dt, but higher derivatives differ). Let f be 
outer normal and r pointing away from the corner [clockwise from v direction). 

The angle between radial and tangential direction is < C{R'Y , because of 
(|5.1.1[) . We can make R' small to control it. Note that C can change from 
line to line. 

w,>{l- C{R'r)r~^w^ - C{R'r\wr\ 

= (1 - CiR'yy-^fiBe-^'l' - C{R'yXr^-\l - fie-^'f') 
= {B^ie-^'f' ^C{R'y)r^~^ 
\wr\ < (1 - C{RT))\wr\ + C{R'rH\ < (A + C{R'rB^i)r^~' 

^^We need full ellipticity later, but here only a'"'', a'f"f' > matter. 
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For any constant T, 

+ Twr > (^Bfie-^'^ - C{R'y - T(A + C{R'y Bfi)y^-^ . (5.1.14) 



D'^u and D{g^) 



{g^\r.= d.,.{g\x,u,Du{x))) = gl,u,k + C (j = l,2). 

(C is a remnant of terms containing only u{x),Du(x) which can be estimated 
by Cu via (|5.1.3p .) Combined with (|5.1.4p we have a system 



9,^ 9,^ 

. ffpl ffp2 





Wll 









Wl2 




(.9^)1 




.""22. 







c. 



The system matrix has determinant 



a^\gl.f-2a^^gl.gl.W\gl.f = ((.g^)^)^A(.gi)^ ™ 



^'^ > (5pC, >0. 



Combined with an upper bound on the matrix (from (|5.1.6p ) we obtain that 

D{x)SI{g^)+C (5.1.15) 



Mil 
U22 



where D{x) G IR^'^ with |i:>(a;)| < C. 
Interior PD inequality 



= g\ka}^Uki] + (ffpfcpi^ijf + C')itifc + Ma^^Wij + C. (5.1.16) 

Again we use C to abbreviate terms depending on u, Du only, which are bounded 
by (l5Xa . 

Take ^fc of (j5T4)l : 

a'-^Utjk + {a^iUki + C)wij = 
Multiply with g^^ and substitute into (|5.1.16p : 

a'^Wy = -(C + glka_^iUM)uij + {glk^iUji + C)Mifc + Ma^^Wij + C. 

Apply (15.1.151) : 

a*^«,, = 7\/a*^u;y + e*^(5i),(5i), + e\g% + C 
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for some coefficients e'^^{x), e*(x) = C. Use Vi = {g^)i + Mwi to obtain 
< r_e--BV-2 + CM2|Vwp + C 



<i 



We take M = M'{R')-^ with A/' > 1 to be determined 



a ■'vij +q Vi < 



ii) 



Given any M' and R' £ (0,1], we take fi (and A, which depends on it) so 
small(er) that the (///) term is dominated by the (/) term. Then we pick A 
even small(er) so that the (//) term is also dominated by the (/) term. In both 
cases we use that {r/R')^ > {r/R'f^ because Q < r < R' . Now 

a'^v.j + q'v, < -Ml^e-^^{r/R')\-^ + C. 

>(i?/)-2 

Finally we take R' so small that the C term is dominated. We obtain 

a'^Vij + q'vi < inU r\ Br'{xo). 
Therefore v cannot have minima in that set. 
Boundary PD inequality 

Take dr of the boundary condition (|5.1.5p on to obtain 

ffpxUrr + gl^U^r + C = 0. 

Wc plan to take a linear combination of this with the interior equation (|5.1.4p 
to obtain a new boundary condition of the form 



= .gp^W^i/ + {Tg^u + glr)Urv + Tg^^Urr 

This problem can be written 



(5.1.17) 



[yi 2/2 -1] 







0. 
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It is solvable if and only if the determinant is 0, which is a linear equation for 
T: solution 

y_ C ^ C 

a""" ■{glrgl.-gl.gl.) a-- AeiG' 

The denominator is bounded away from by (|5.1.7p and (|5.1.9p . so |r| < C. 
The solutions are 

^ ff^ ^ a^^{gl.f-2a^-gl.gl.+a--{gl.)\ 

a""^' a^-'detG 

clearly |yi|, \y2\ < C. Therefore we are justified in writing C on the left-hand 
side of (15.1.171) . 

Now: 

V, + Tvr = M{w, + Twr) + {g^)u + T{g% 

MifiBe-^'l' - C{R'y - T(A + CiR')" B^l))r^-^ + C 



15.1.171 

= M'{R')-^(fiBe-^^ - C{R'Y - T(A + C{R'Y B^i))r^-^ + C 



We can choose A and R' so small that 

V, + Tvr > M'{R')-^^r^-^e-^'f' + C for all r € (0, R'). 
The C term is bounded, so we can pick R' so small(er) that 
+Tvr>0 for all r £ (0, R'). 

This inequality implies that v cannot have minima on H Bn' (x) because 
otherwise < and Vr = in the minimum point would cause a contradiction. 

Conclusion of corner C^'^ 

Let 

/' := : </) € U' := {(r,0) : r G {0,R'), <t> S [r) , <t>^ [r))} . 

We are still free to choose Ad' (> 1). Since Mw > ^ on /', whereas Du{x)) = 

C, we can take M' so large that u > on /'. 

V must attain its minimum over U' on /' U {.to} U F^. On the boundary 
condition g^ = yields v = w > 0. Therefore v attains its minimum u = in 
xq. Hence v — g^ + Mw > in J7, so 

y - - 2 \R' ' 
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Moreover we can apply the same arguments to v := — Mw, by switching 
inequahties and signs and replacing "minimum" with "maximum" whereever 
appropriate. Then 




The same arguments, but with the roles of and reversed, yield a bound 
on g^. Now we argue that g^ and 5^, as functions of Vu with x and u{x) held 
fixed, are uniformly functionally independent, by (|5.1.9p . Therefore 

|Vu(a;) - Vu(a;o)| < Cr^. 



Pull C^'" result 

This concludes the proof of (|5.1.10p . The C^^_^ result follows from |Lie88a[ 
Corollary 1.4], where R is the distance from the corner, and osco^, Du/R = i?'^^^ 
by (|5.1.10p . □ 



5.2 Free boundary transformation 

For a column vector x and a column- vector- valued function u, let be the 
transpose of the gradient and u^x the Hessian. For "row" instead of "column" 
let Ux be the gradient. Boundary normals are row vectors. 

Proposition 5.2.1. Consider a coordinate transformation y = y{x,u{x)). 

1. It is nondegenerate if and only if 

Vx + VuUx 

is regular. Let Cd he an upper bound for the norms of the matrix and its 
inverse (for some fixed matrix norm). 

2. A quasilinear equation 

A{x, u{x),Ux{x)) : Uxx{x) + b{x, u{x), Ux{x)) = 

transforms into another quasilinear equation. The first equation is elliptic 
if and only if the second one is, and the ellipticity constants are compara- 
ble, up to constant factors depending only and continuously on Cd . 

3. Consider a boundary condition 

g{x,u{x),Ux{x)) = 0. 
As before we say it is oblique, with obliqueness constant 60 > 0, if 

\ngp\ > So\gp\, 
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where p represents Ux- If 

D:=\l- Uxivx + yuu^y^yul > So > 0, (5.2.1) 

then the new boundary condition is oblique if and only if the old one is, and 
the obliqueness constants are comparable, up to constant factors depending 
only and continuously on So and Cd- 

4- Under the same assumptions, consider two functions g^,g^ of Ux- The 
smaller angle between g^ , g^ is nonzero if and only if the smaller angle 
between g^^g^ is nonzero. They are comparable, up to constant factors 
depending only and continuously on Cd,S]j. Here q represents Uy. 

Proof. 1. Obvious. 

2. The new top-order coefRcient matrix is 

iVx + yuUx)A{yx + yuUx)^ 

It is still symmetric positive definite, and the cUipticity constant bound is 
obvious. 

3. We omit the x and u{x) arguments from the notation of the boundary 
conditions. 

u{x) = u{y{x,u{x))). 

Take the x derivative: 

Ux = UyijJx + yuUx) ^ Uy = Ux{yx + yuUx)~^ ■ 

If h{uy) = is the boundary condition in y coordinates, then 

giUx) = h{ux{yx + VuUx)~^)- 

Expressed with p, q: 

dip) = h{p{yx +yuP)~^), 

so 

gp = [ivx+Vupy^-pi.yx+VupY^yui.yx+yupy^^hq = {i-p{yx+yupy^yu){yx+yup)~^hq. 

With p = Ux 'we have 

gp (1 - Ux{yx + yuUx)~^yu){yx + yuUx)~^hq. 

Let n be the normal in x coordinates, N the normal in y coordinates. 
Then 

N{yx + yuUx) 

n — - — --. 

\N{yx + yuUx)\ 
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So 

nOr, = : — ; r; N n„. 

|iV(y, + y„u,)| ' 

Clearly ng^ 7^ iff Nhq ^ 0, as long as Z) > 0. The constant comparison 
is obvious. 

4. The two functions transform like the boundary condition above: 

5p = (1 - u^{yx + yuUx)~^yu){yx + yuUx)~^hg 

where g''{ux) = h^{uy). Then clearly and g^ are collinear if and only 
if /ij and h?^ are, and the constants are obvious. 

□ 

Remark 5.2.2. We use the Sherman-Morrison formula: for an invertible matrix 
M, the rank 1 perturbation M + ah^ is invertible if and only if b'^M~^a ^ —1. 
Taking M = yx + yuUx, a = —yu and b — , we see that (|5.2.ip is precisely 
the condition for to be invertible as well. 
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